Homogenization of Elliptic Boundary Value Problems 

in Lipschitz Domains 

Carlos E. Kenig* Zhongwei Shen^ 

Dedicated to the Memory of Bjorn Dahlberg 

Abstract 

In this paper we study the LP boundary value problems for C(u) = in R^j" 1 , where 
C = — div(v4V) is a second order elliptic operator with real and symmetric coefficients. 
Assume that A is periodic in x^+i and satisfies some minimal smoothness condition in 
the Xd+i variable, we show that the L p Neumann and regularity problems are uniquely 
solvable for 1 < p < 2 + 5. We also present a new proof of Dahlberg's theorem on the 
LP Dirichlet problem for 2 — 5 < p < oo (Dahlberg's original unpublished proof is given 
in the Appendix). As the periodic and smoothness conditions are imposed only on 
the Xd+i variable, these results extend directly from to regions above Lipschitz 

graphs. Consequently, by localization techniques, we obtain uniform LP estimates for 
the Dirichlet, Neumann and regularity problems on bounded Lipschitz domains for a 
family of second order elliptic operators arising in the theory of homogenization. The 
results on the Neumann and regularity problems are new even for smooth domains. 

MSC 2000: 35J25. 

Keywords: Boundary value problem; Periodic coefficient;Lipschitz domain. 



1 Introduction 

Let C = — div(AV) be a second order elliptic operator denned in ~§t d+1 = {X = (x,t) E 
M. d x R}, d > 2. We will always assume that the (d+ 1) x (d+ 1) coefficient matrix 

A = A(X) = (ajj(X)) is real and symmetric, (1-1) 

and satisfies the ellipticity condition, 

/^| 2 <MX)£4-<-|£| 2 forallX,£GM d+1 , (1.2) 
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where \x > 0. In this paper we shall be interested in boundary value problems for C(u) = 
in the upper half-space = R d x (0, oo) with LP boundary data, under the assumption 
that the coefficients are periodic in the t variable, 

A(x, t + 1) = A(x, t) for (x, t) £ lR d+1 . (1.3) 

More precisely, we will study the solvabilities of the LP Dirichlet problem (D) p 

C(u) = inn = Hf 1 , 

u = f £ L p (dVt) n.t. on dVL and (it)* £ L p (dVt), 1 ' } 

the LP Neumann problem (N) p 

C{u) = in tt = M. d +\ 

dv (1-5) 
_ = / e L p (dn) on <9fi and iV(Vw) £ L p (<9ft), 

where |^ denote the conormal derivative associated with operator £, and the L p regularity 
problem (R) p 

' C(u) = in n = R d +\ 

• 1 (1-6) 
u = / e W^ p {dn) n.t. on 9fi and N(Vu) G L p (9fi). 

Here («)* denotes the usual nontangential maximal function of m and iV(Vw) a generalized 
nontangantial maximal function of Vw. By u = f n.t. on <9fi we mean that u(X) converges 
to f(P) as X —>■ P nontangentially for a.e. P £ Under the periodic condition f 1 1.3 1) as 
well as some (necessary) local solvability conditions on C, we will show that the LP Dirichlet 
problem is uniquely solvable for 2 — 5 < p < oo, and the LP Neumann and regularity 
problems are uniquely solvable for 1 < p < 2 + 5. Furthermore, the solution to the Dirichlet 
problem satisfies the estimate ||(w)*llp — C|l M llp5 while the solutions to the LP Neumann and 
regularity problems satisfy ||-/V(Vm)|| p < C||f^||p and ||iV(Vit)||p < C||V x m|| p , respectively. 
These results extend the work on the LP boundary value problems in M^ +1 for second order 
elliptic operators with t-independent coefficients (or in bounded star-like Lipschitz domains 
for operators with radially independent coefficients). As we shall discuss later, since the local 
solvability conditions may be deduced from certain minimal smoothness conditions in the t- 
variable and our periodic condition is only imposed on the t-variable, the LP global estimates 
on (it)* and iV(Vu) extend directly from M^ 1 to the regions above Lipschitz graphs. As a 
consequence, by well known localization techniques, we obtain uniform estimates for the LP 
boundary value problems in bounded Lipschitz domains for a family of second order elliptic 
operators arising in the theory of homogenization. 

We should point out that the result mentioned above for the Dirichlet problem is in fact 
due to the late B. Dahlberg [9] (unpublished, personal communication). His proof, which 
is given in the Appendix for the sake of reference, depends on the ingenious use of Green's 
functions and harmonic measures. As in the case of his celebrated theorem on the Dirichlet 
problem for Ait = in Lipschitz domains [TQl [TTJ , Dahlberg's method does not extend to the 
Neumann and regularity problems. Motivated by the work of Jerison and Kenig [T6 t [T71 [T8] . 
we then seek to establish the global L? Rellich type estimate, 

ii du.. ,._ h „ x 

bH 2~ V S U 2 (1-7) 

ov 
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for suitable solutions of elliptic operators with t-periodic coefficients. In Section 3 we develop 
some new integral identities which play a crucial role in the proof of (11.71) and which may 
be regarded as the Rellich identities for operators with t-periodic coefficients. Indeed, in 
the case of constant or t-independent coefficients, Rellich identities are usually derived by 
using integration by parts on a form involving The basic insight here is to replace the t 
derivative of u by the difference Q(u)(x,t) = u(x,t + l) —u(x,t). The t-periodicity of A(x,t) 
is used in the fact that Q commutes with C. In particular, Q(u) is a solution whenever 
u is a solution. We further remark that these integral identities allow us to control the 
near-boundary integral 



by the integral of |§^| 2 or (V^w] 2 on M. d . The desired estimate (11.71) then follows by local 
solvability conditions. 

With (11.71) at our disposal, we are able to give another proof of Dahlberg's theorem on 
the Dirichlet problem in Section 4 and more importantly, solve the L 2 regularity problem in 
Section 5 and the L 2 Neumann problem in Section 8. Furthermore, using the L 2 estimates 
and following the approach developed in [13] and [21], we establish the solvability of the LP 
regularity and Neumann problems for 1 < p < 2 in Sections 6 and 9. Although the range 
2 < p < 2 + 5 may be also treated by the real variable arguments used in [131 EI] , we choose 
to use a relatively new real variable approach, based on the weak reverse Holder inequalities 
[251 [26] . I n particular it allows us to show that for elliptic operators with coefficients satisfying 
(ll.ip - (ll.2l) . if the L 2 regularity problem (R)2 is solvable and ^ + - = 1, then the Dirichlet 
problem (D) g and the regularity problem (R) p are equivalent. 

There exists an extensive literature on boundary value problems with minimal smoothness 
assumptions on the coefficients or on the boundary of the domain in question. In particular 
the L p Dirichlet, regularity and Neumann problems for Au = in Lipschitz domains, which 
are closely related to the problems we investigate here, were well understood thanks to 
[TOl [TT1 fl~6l [T7l [T8l l29l fl~3] . Deep results on the L p Dirichlet problem for a general second order 
elliptic equation C(u) = may be found in [HI [El [HI EE] (see [19] for further references). 
The LP regularity and Neumann problems for C(u) = were formulated and studied in 
[21~1 [22]. If the coefficient matrix A is t- independent, the L p solvability of the Neumann 
and regularity problems in the upper half-space was essentially established in [21] for the 
sharp range 1 < p < 2 + 5, although the paper only treats the case of radially independent 
coefficients in the unit ball. We also mention that the L 2 boundary value problems for some 
operators with t-independent complex coefficients was recently studied by the method of 
layer potentials in [I]. Related work on L p boundary value problems for operators with 
t-independent, but non-symmetric coefficients may be found in [20l [23] . Regarding the 
conditions on the t variable, it worths mentioning that the global estimate ||(w)*|| p < C||/||p 
may fail for any p, even if the matrix A(x,t) satisfies fll.ip - fll.2D and is C°° in M. d+1 . 

Since the local solvability condition on £ may be deduced from certain minimal smooth- 
ness assumption in the t variable on the coefficient matrix A, as we alluded to earlier, 
Dahlberg's theorem on the Dirichlet problem and our results on the Neumann and regular- 
ity problems extend readily to the case where 




(1.8) 



D = {{x,t) e R d+1 :t>i){x 



)} 



(1.9) 
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is the region above a Lipschitz graph, using the bi-Lipschitzian map (x,t) — > (x, t — ip{x)) 
from D to Indeed, let 

7]{p) = sup {\A(x, si) - A(x, s 2 )| : x £ IR d and |s x - s 2 | < p) (1.10) 

and assume that 



1 Hp)) 



2 



■dp < oo. (1-11) 

/o P 

The following three theorems are proved in Section 10. The constant 5 in Theorems 11.11 E2] 
and 11.31 depends only on d, p, rj(p) and || V^lloo- 

Theorem 1.1. Let £ = — div(A'V) with coefficient matrix satisfying ( ti.il) . ( ti.jjj) . ( ti.gj) and 
M.ll\) . Let D be given by U.9\) . Then there exists 5 £ (0, 1) such that for any f £ L p (dD) 
with 2 — 5 < p < oo, there exists a unique solution to C(u) = in D with the property 
that (u)* £ L p (dD) and u = f n.t. on dD. Moreover, the solution u satisfies the estimate 
\\{u)*\\ p < Cp\\f\\p, where C p depends only on d, p, p, p(p) and || V^Hoo- 

Theorem 1.2. Under the same assumptions on C and D as in Theorem \l.l\ there exists 
5 > such that given any g £ L p (dD) with 1 < p < 2 + 5, there exists a solution, unique up 
to constants, to C(u) = in D such that N(Vu) £ L p (dD) and |^ = g on 3D. Moreover, 
the solution satisfies the estimate ||A^(Vm)|| p < Cp||g|| p; where C p depends only on d, p, p, 
r](p) and HVxVlloo- 

Let VtanU denote the tangential derivatives of u on dVL. We say / £ W 1,p (dVt) if / £ 
Lf oc (dn) and V tan f EL p (dn). 

Theorem 1.3. Under the same assumptions on £ and D as in Theorem \l.l\ there exists 
5 > such that given any f £ W 1,p (dD) with 1 < p < 2 + 5, there exists a unique solution 
to C{u) = in D such that N(Vu) £ L p (dD) and u = f n.t. on dD. Moreover, the solution 
satisfies the estimate \\N(Vu)\\ p < C P WVtanf\\ P , where C p depends only on d, p, p, p(p) and 

HV^IIoo- 

As we indicated earlier, Theorem II .11 as well as Theorem 1 1 . 41 b elow are due to B. Dahlberg 
[9]. However, under the i-periodic condition, Theorems ll.2l and ll.3l are new even for operators 
with smooth coefficients on smooth domains. We further point out that by well known 
localization techniques, Theorems 11.21 and 11.31 allow one to establish uniform estimates 
for the LP boundary value problems in Lipschitz domains for a family of second order elliptic 
operators arising in the theory of homogenization (see e.g. [7]). 

More precisely, we consider 

C £ = -div [A (—) V ) , e > 0, (1.12) 



e J 

where the matrix A(X) is periodic with respect to the standard lattice, 

A(X + Z) = A(X) for any X £ R d+ \ Z £ Z d+1 . (1.13) 
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In view of (11. lip , we will assume that A(X) is continuous and its modulus of continuity 
satisfies 

{UJ{P))2 dp < oo, (1.14) 



lo P 

where u(p) = sup {\A(X) - A{Y)\ : X,Y eR d+1 and \X-Y\<p). 
The proofs of the following three theorems are given in Section 11. 

Theorem 1.4. Assume that the coefficient matrix A(X) is real, symmetric, and satisfies 
the ellipticity condition U.S\) . the periodicity condition U.13\) and the smoothness condition 
(fTTT^p. Let Q be a bounded Lipschitz domain in M. d+1 . Then there exist constants 5 G (0, 1) 
and C > independent of e, such that if f G L p (dQ) for some 2 — 5<p< oo, then the unique 
solution to the LP Dirichlet problem: C £ (u £ ) = in O, u £ = f n.t. on dQ, (u e )* G L p (dQ), 
satisfies the estimate ||(w e )*|| p < (7||/|| p . 

Let J^ 1 denote the conomal derivative associated with the operator C £ . 

Theorem 1.5. Let Q be a bounded Lipschitz domain in M. d+1 . Assume that A(X) satisfies 
the same assumption as in Theorem \l-4\ Then there exist constants 5 G (0,1) and C > 
independent of e, such that given any g G L p (dQ) with 1 < p < 2+5 and mean value zero, the 
unique solution to the L p Neumann problem: C £ {u £ ) = in f2, |^ = g on <9f2, N{SJu £ ) G 
L p (dfl), satisfies the estimate \\N(Vu £ )\\ p < C\\g\\ p . 

Theorem 1.6. Let Q be a bounded Lipschitz domain in M. d+1 . Assume that A(X) satisfies 



the same assumption as in Theorem L4_. Then there exist constants 5 G (0, 1) and C > 
independent of e > 0, such that given any f G W 1,p (dQ) with 1 < p < 2 + 5, the unique 
solution to the L p regularity problem: C £ {u £ ) = in Q, u £ = f n.t. on dQ, N(Vu £ ) G 
L p (dQ) satisfies the estimate \\N(Vu e )\\ p < C { || V tan /|| p + \d£l\^ \\f\\ P }. 

We remark that Dahlberg's work on operators with periodic coefficients was inspired 
by a series of remarkable papers [21 El IH El [6] of M. Avellaneda and F. Lin on elliptic 
homogenization problems. In particular, the uniform estimate ||(w e )*|| p < i n Theorem 

11.41 for 1 < p < oo was established on C 1,7 domains by Avellaneda and Lin in [21 El E] for 
second order elliptic equations and systems in divergence form with C a periodic coefficients. 
Related work on the boundedness of Riesz transforms V(£ e ) -1 ^ 2 on Lipschitz and C 1 domains 
may be found in [2H] . However, to the authors' best knowledge, the uniform LP estimates for 
periodic operators have not been studied before for the regularity and Neumann problems; 
Theorems 11.51 and 11.61 are new even for smooth domains. 

Finally we mention that as in the case of operators with constant coefficients, the Rellich 
estimates we develop in this paper, can be used to solve the L p boundary value problems by 
the method of layer potentials. This would enable us to represent the solutions in Theorems 
ll.4[ 11.51 and 11.61 in terms of layer potentials with density functions uniformly bounded in 
appropriate spaces. Even more importantly, the method of layer potentials can be applied 
to elliptic systems and would allow us to establish uniform LP estimates for elliptic systems 
with periodic coefficients in nonsmooth domains. This line of ideas will be fully developed 
in a forthcoming paper [24J. 
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2 Notations and Preliminaries 

We will use X, Y, Z to denote points in M. d+1 and x, y, z points in M. d . Balls in ~R d are denoted 
by B(x,r) = {y G M. d : \y — x\ < r}. By T(x,r) we mean the cylinder B(x,r) x (0, r). We 
will use V for the gradient in M. d+1 and V x for the gradient in M. d . 

The classical Dirichlet problem in Q = . 

Let / G C(R d ) such that f(x) — * as \x\ — > oo. Then there exists a unique weak solution 
to £{u) = in R^. +1 such that u G C(R d+1 x [0,oo)), u(x, 0) = f(x) for x G R d and u is 
bounded in R^ +1 . Moreover, the solution satisfies sup{|w(x,t)| : x G R d ,t > 0} < ||/||oo- 
The uniqueness follows from Louiville's Theorem by a reflection argument. To establish the 
existence, one may assume that / > 0. Let Qk = T(0, k) = -8(0, k) x (0, k) and (p a continuous 
decreasing function on [0, oo) such that < (p < 1, (p(s) = 1 for < s < (1/2) and <p(s) = 
for s > (3/4). Let Uk be the solution to the classical Dirichlet problem for C(u) = in 
T(0,k) with boundary data f(x)if(\x\/k) on {(x, 0) : \x\ < k} and zero otherwise. Since 
Uk < Mfe+i on Qk and ||wfc||L°°(n fc ) < ll/lloo, it follows that u(x) = limfc^oo Uk{x) exists and is 
finite for any x G IR+ +1 and llwll^^d+i) < ||/||oo- Also, we may deduce from the estimate 

sup \uk — u m \ < C{uk — u m )(0, £/2), fork>m>£ 
n t 

that as k — > oo, Uk converges to u uniformly on Qi for any i > 1. This implies that C(u) = 
in n = u is continuous in ft, and u(x, 0) = f(x). 

By the Riesz representation theorem, there exists a family of regular Borel measures 
{u x : X G on R d , called £-harmonic measures, such that 

u(X)= [ f(y)du; x (y). (2.1) 

Let Vk be the solution to the classical Dirichlet problem in with data cp(\x\/k). It 

follows from the uniqueness that Vk(X) — > 1 as A; — > oo for any X G This shows 

that uj x is a probability measure, as in the case of bounded domains. Furthermore, let uj x 
denote the ^-harmonic measure for the cylinder -8(0, k) x (0, k). It is not hard to see that 
uj x (B) -> lu x (B) as k -> oo, for any ball Bcl d and X G E^ +1 . 

Green's function on fi = 

Let Gk{X,Y) denote the Green's function for £ on the domain Qk = T(0,k). Since 
< G k (X, Y) < G k+ i(X, Y) for X, Y G Q k , the Green's function for £ on = may be 
defined by 

G(X, Y) = lim G fe (X, Y) for X, F G R+ +1 and X ^ Y. 

k^oo 
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By well known properties of G^, we have 

C Ct a 
0<G(X,Y)< lx _ Yri and G(X, Y) < ^ _ - {d _ 1+a (2.2) 

for X = (x, t), Y G R+ + \ where C and a are positive constants depending only on d and \x. 
Since u£(B(y,r)) « r d " 1 G fc (X, (j/,r)) for any X G fl*. \ T(y,2r) [8J, by letting fc — > oo, we 
obtain 

u x (B{y, r)) « r^GpT, (y, r)) for any X £ 1% 2r). (2.3) 

Two properties of nonnegative weak solutions. 

Let u be a nonnegative weak solution of C(u) = in T(xq, 2r) for some Xq G M. d and 
r > 0. Suppose that u G C(T(xo,2r)) and u(x,0) = on B(xo,2r). Then we have the 
boundary Harnack inequality, 

u(x, t) < Cu(x , r) for any (x,t) G T(x ,r). (2.4) 

We will also need the following comparison principle, 

— — < C— for any (x,t) G T(x ,r), (2.5) 

t>(x, r) f(xo,r) 



where u, v are two nonnegative weak solutions in T(x ,2r) such that u, v G C(T(xo,2r)) 
and 0) = t>(x, 0) = on -B(xo, 2r) (see e.g. [8]). The constants in ( 12. 4ft and ( 12.51) depend 
only on (i and //. 

Difference operator Q. 

Define 

Q(«)(x,t) =u(x,t + 1) -u(x,t) for (x,t) ER d+1 . (2.6) 
Clearly, Q(J^) = ~£rQ( u ) f° r i = 1, • • • , d + 1. It is also easy to verify that 

QM-QWQM = «Q(<;) + <?M (2.7) 

and 

Q(u)(y, s) dyds = / / u(y,s)dyds— / / u(y,s)dyds. (2.8) 

T(z,t) it J B{x,t) JO JB(x,t) 

Nontangential maximal functions. 

For a function u defined in we let 

3t „ \ 1/2 

2 



jV(u)(z) = sup <( ( ^ ^ jf^ ^ \u(y,s)\ 2 dydsj : < t < oo } . (2.9) 

This is a variant of the usual nontangential maximal function (u)* which is defined by 

(u)*(x) = sup {\u(y,t)\ : < t < oo and \y - x\ < t}. (2.10) 
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Clearly, N(u)(x) < (u)*(x). If u is a weak solution to C(u) = in R+ +1 or any function that 
has the property 



\u(x,t)\ < C ( —j-r / / \u{y,s)\ 2 dyds 

\r d+1 JJB(( x ,t),r) 



1/2 



whenever B((x,t),r) C then ||A^(m)|| p « ||(w)*|| p for any < p < oo. We also need to 

introduce N r (u) and (it)* which are defined by restricting the variable t in (12.91) and (12.101) 
respectively to < t < r/2. The definitions of N(u) and (u)* extend naturally to regions 
above Lipschitz graphs and to bounded Lipschitz domains. 

We conclude this section with a lemma which allows us to estimate \\N(Vu) \\ p by || A^(Vu) || p 
and \\{Q(u))%. 



Lemma 2.1. Let u G C 1 (T(0, 3R)) be a solution to C(u) = in T(0, 3R) for some R > 10. 
Let f(x) = u(x, 0). Then for any x G B(0, R), 

N R (Vu)(x) <Cm([N4Vu) + |V X /| + {Q(u)T r ]xb(o,3R)){x), (2.11) 

where M denotes the Hardy-Littlewood maximal operator on M, d and C depends only on d 
and \i. 

Proof. We begin by fixing (x , t ) G R d+1 with x G 5(0, R) and t Q G (2, R/2). Note that by 
Cacciopoli's inequality and interior regularity of weak solutions, 



1 



£ J\s-to\<% JB(x ,t ) 



<t(t^I I Hy,s)-E\*dydsY (2-12) 

J \t Q J\s-t \<aito JB(x ,Pito) J 

< / / \u{y, s) - E\ dyds, 

t J\s-t \<a 2 t J B(x ,j3 2 t ) 

where (1/2) < a x < a 2 < (3/4), 1< /3 X < (3 2 < (3/2) and £ G K. To estimate the right-hand 
side of (12.121) . we use 

\u(y, s)-E\< \u(y, s) - u(y, 0)| + \u(y, 0) - E\ 
and choose E to be the average of u(y, 0) over B(x , Pit ). By Poincare's inequality, we have 

i / / \u(y, 0) - E\dyds 

<S / iV^/ldx^CAfdVx/lxB^))^). 

J JB(xo,/3 2 to) 

Next we observe that for any s G (0, 2t ), 

K2/ ) s)-w(y,0)|<Ct (Q(w))i(|/)+ / |Vu(y,t)|dt. (2.14) 



■ 1/2 

| Vu{y, s)\ 2 dyds 



1/2 



In view of ( 12.121) . the first term in the right-hand side of (12.141) can be handled easily by 
M ^(Q(u))* r xb(o,3R)) (%o)- Thus it remains to estimate 



£q JB{x ,j3 2 t ) Jo 



1 

\Vu(y, t)\ dtdy. 



This may be done by observing that 



C [ N 4 (Vu)dx> f f \Vu(y,s)\dyds. 

JB(x Q ,4t ) JO JB(x ,2t ) 



□ 



Remark 2.2. It follows from (12.111) that for 1 < p < oo, 

/ \N R (Vu)\ p dx<C [ \V x u{x^)\ p dx + C [ \N 4 (Vu)\ p dx 

JB(0,R) JB(0,3R) JB(0,3R) 

+ C f \(Q(u)) R \ p dx 

JB(0,3R) 

for any R > 10. Thus, if u G C^R* x [0, oo)) and C{u) = in we may let R -> oo in 

fl2TT5|) to obtain 



(2.15) 



\N(Vu)\ p dx < C I \S7 x u(x^)\ p dx + C \N A (Vu)\ p dx 

jRd (2.16) 
+ C \(Q(u))*\ p dx, 



where C depends only on d, p and \x. 



3 Integral identities 

In this section we develop some new integral identities that will play a key role in our 
approach to the L 2 boundary value problems in R+ +1 for elliptic equations with t-periodic 
coefficients. 

Lemma 3.1. Let £ = —div(AV) with coefficients satisfying conditions / tl.lj) -/ TOj) . Let 
Q C R d+1 be a bounded Lipschitz domain and u G W 1,2 (Vt) a variational solution to the 
Neumann problem: C(u) = in O, |^ = g G L 2 {dVl). Assume that u G W 1,2 {Vli) where Vl\ 
contains Q and {(x, t + 1) : (x, t) G Q}. Then 

Jn ux i ux j Jn ux t ux j ^2 j\ 

= —2 g Q(u) da, 
Jan 

where indices i,j are summed from 1 to d+ 1 and we identify Xd+i with t. 
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Proof. Since Q(u) £ W 1,2 (£l), it follows from the definition of variational solutions that 

f f du du 

/ gQ(u)da= / a^— ■ Q(—) dX 

Jon Jn c>x j ® x i 

2 In ° lJ ' { 9xj ^ ^ dxi ^ dxi ® ^ dxj ^ } 

where we have used the symmetry condition = aji in the second equation. In view of 
(12. 7p . this gives 

from which the desired identity follows. □ 

In the rest of this section we will assume that satisfy the periodic condition (jl .3j) . To 
justify the use of integration by parts, we will also assume that both eiy and u are sufficiently 
smooth. However all constants C in this section depend only on d and \x. 

The following lemma is crucial in our argument for the L 2 Neumann problem. 



Lemma 3.2. Let £ = — div(AV) with C°° coefficients satisfying U.l\) . (CQ|) and U.3\) . 
Let R > 10 and u G C 1 (f2) be a solution to C(u) — in O = T(xq,3R) and |^ = g on 
B(x ,3R) x {0}. Then 

2 . C 



Vu\ 2 dxdt<C / \g\ dx + — \Vu\ 2 dxdt, (3.2) 

'0 JB(x ,R) Jb(x 0i 2R) R J Jt(x ,3R)\T(x ,R) 

where C depends only on d and \i. 

Proof. We may assume that x = 0. Let f2 r = T(0, r) = B(0, r) x (0, r) for r 6 (_R, 2_R). By 
the periodicity assumption (11.31) and H2. 81) . we have 

r r du du f f du du 

- JL a " Q{ ^ ■ ar? dxdi = - JL Q[a "^ ' ar? dzdt 

du du , , f r+1 f du du , , 

claa— — • — — dxdt — I an— — ■ — — dxdt. 



'0 J B(0,r) J r JB(0,r) &Xj 

It then follows from (13. ip that 

f f du du f f du du 

/ / aij—-—dxdt+ // a ij Q(—)Q(——)dxdt 

JO J B(0,r) OXi OXj J JQ r OXi OXj 

=- 2 L^ Q{u)d ° + l 



r+1 f du du 



ldQ r UL/ Jr JB(0,r) 3 dXi dxj 



(3.3) 
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for r G (R, 2R). By the Cauchy inequality, this leads to 

/ / \Vu\ 2 dxdt < 5 f \Q(u){x,0)\ 2 dx + C s [ \g(x)\ 2 dx 

J0 JB(0,r) JB(0,r) J B(0,r) 

+ C \Vu(x,r)\ 2 dx + C / \Q(u)(x,r)\ 2 dx 

JB(0,r) JB(0,r) 



B(0,r) JB(0,r) 

+ C I I \Vu\ 2 dxdt + C I I \Vu\ 2 dadt 

B(0,r) JO JdB(0,r) 



(3.4) 



C I / \Q(u)\ 2 dadt, 

'0 JdB(0,r) 



for any < 8 < 1. Using 



|Q(w)(a;,t)| 2 rfa; < / / |Vu| 2 dxds 

B(0,r) Jt J B(0,r) 

rr+1 f 

2 j_Jj. ^ \ I IV7„.|2 



\Q(u)\ 2 dadt < / / |VM| 2 do-(it, 

JdB(0,r) Jo JdB{0,r) 



we obtain from (13.41) that 



[ \Vu\ 2 dxdt<C [ \g{x)\ 2 dx + C [ \Vu{x, r)\ 2 dx 

Jb(0,R) Jb(0,2R) Jb(0,2R) 

rr+1 /> 

C / \Vu\ 2 dxdt (3.5) 

Jr Jb(0,2R) 
3R 

C l l \Vu\ 2 dadt, 



JdB(0,r) 

for r G (R,2R). The desired estimate now follows by integrating both sides of (13. 5p with 
respect to r over the interval (R, 2R). □ 



Remark 3.3. Under the same assumption as in Lemma I3.2[ we also have 



2 . C 



\Vu\ 2 dxdt<C / \g\ 2 dx + — \Vu\ 2 dxdt. (3.6) 

JB(x ,R) Jb(x ,2R) R J Jt(x ,3R)\T(x ,R) 

To see (13.61) . one simply replaces Q{u) with Q{u) = u(x,t + 6) — u(x,t) in the proof. 

The following lemma is the key in our approach to the L 2 regularity problem. 

Lemma 3.4. Under the same assumptions on C as in Lemma UT^, there exists C = C(d, n) > 
such that 

[ [ \Vu\ 2 dxdt<C [ \V x f\ 2 dx+ ^ I! \Vu\ 2 dxdt, (3.7) 

Jo Jb(x ,R) Jb(x ,2R) R J Jt(x ,3R)\T(x ,R) 



where u G C 1 (T(xq, 3R)) H C 2 (T(x , 3-R)) is a solution to C{u) = in T(xq, 3R) and u = f 
on B(x ,3R) x {0}, R > 10. 
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Proof. We may assume that Xq = 0. As in the proof of Lemma l3~2l we let Q r = B(0, r) x (0, r) 
for r G (R, 2R). By the periodicity assumption, Q(u) is a solution in Q r . It follows that for 
r G (R,2R), 



/" 9 

/ u ■ —Q(u) da 
Jdu r dv 



dflr 

un i' f -zr \oij(x,t + s)-^-(x,t + sH dsda(x,t) 
Jo 9s { dxj 

urii-l^- [ aa(x,t + s)^-(x,t + s) ds da(x,t) 
[dtJo dx-j 1 



an r 

u 



s )tt~~( x > t + s )ds )■ da(x, t). 

OXj 



(3.8) 



d d f f 1 du \ 

'y li Qt~ nd + 1 Q^J jy aij(x,t + s) — (x,t + s)dsj da(x,t) 

Consequently, by the Cauchy inequality, we obtain 

I f |^ • Q(u)da\ 
1 Jan r dv 1 

<C / / \\/tanu(x,t)\ \\/u(x,t + s)\ da(x 1 t)ds 

JO JdQ r 

<5 [ [ \Vu(x,t)\ 2 dxdt + C 5 [ \V x u(x,0)\ 2 dx (3.9) 

Jo Jb{0,t) J B(0,r) 

+ C [ \Vu(x,r)\ 2 dx + C [ [ \Vu(x,t)\ 2 dxdt 

J B{0,r) Jr Jb{0,t) 

+ C [ I \Vu{x,t)\ 2 da{x)dt, 

Jo JdB(0,r) 

for any < 5 < 1, where V ta n denotes the tangential gradient on dVL r . This, together with 
f)3.3p . implies that 



n\Vu\ 2 dxdt 
J(0,R) 

<C f \V x f\ 2 dx + C [ \Vu(x,r)\ 2 dx (3.10) 

Jb(0,2R) Jb(0,2R) 



' B(0,2R) J B(0,2R) 

i-3R r rr+i 

12 j n , n \ I iv7„.i2 



+ C / \S7u\ 2 dadt + C / / \Vu\ 2 dxdt, 

Jo JdB(0,r) Jr J B(0,2R) 

for r G (R,2R). The estimate (13. 7j) follows by integrating both sides of (13.101) with respect 
to r over the interval (R, 2R). □ 
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Remark 3.5. By replacing Q(u) with Q(u)(x, t) = u(x,t+6)—u(x,t) in the proof of Lemma 
I3.4[ we may obtain 

I [ \Vu\ 2 dxdt<C [ \V x f\ 2 dx + ^ II \Vu\ 2 dxdt. (3.11) 

JO JB(x ,R) JB(x ,2R) R J JT(x ,3R)\T(x ,R) 

4 The L p Dirichlet problem: a theorem of Dahlberg 

To solve the L p Dirichlet problem, we will impose a local solvability condition on C: for any 
x G R d and < r < 1, 



S B\ 



limsup f^A\ dx<^ ff \u(x,t)\ 2 dxdt, (4.1) 

B(x ,r) t->0 \ t J r J JT(x ,2r) 

whenever u G W 1,2 (T(xo, 4r)) is a nonnegative weak solution of C(u) = in T(xo,4r) such 
that u G C(T(x Ar)) and u(x, 0) = on B(x , 4r). Recall that for Z G \ T(x, 2t), 

^M) „ G((x,t),Z) 
\B(x,t)\ t 1 ' ' 

where uj z denotes the £-harmonic measure on M.+ and G(X, Z) the Green's function on 
R d +\ By applying the condition (JUJ to u{X) = G(X,Z) for Z G \ T(x ,4r), we 

may deduce that uj z is absolutely continuous with respect to the Lebesgue measure on R d . 
Furthermore, the kernel function 

r,/ rrS , ^(5(x,t)) 

if(x, Z) := lim , ^/ 
satisfies the reverse Holder inequality 

1 /2 

' f \K(y,Z)\ 2 dy] <£- [ K(y,Z)dy, (4.3) 



I- 8 ! is / I- 8 ! </B 

for B = B(x,r) with < r < 1, if Z ^ T(x, 2r). In fact, by the comparison principle ( 12. 51) . 
conditions ( 14. ip and ( 14.31) are equivalent. 

The goal of this section is to show that with the additional periodic condition ( II. 3p . the 
reverse Holder inequality (14.31) holds without the restriction r < 1. As a consequence, the LP 
Dirichlet problem for C(u) = in R* +1 is solvable for 2 — 5 < p < oo, under the assumptions 
( Oil . ffi~2l) . fli~3l) and (OI The following theorem is due to B. Dahlberg [9]. 

Theorem 4.1. Let C = — div(AV) be an elliptic operator satisfying conditions U.l\) . U.S\) . 
/ li.gj) and ( |^.i| ). T/ien i/iere exists 5 = 5(d, /i, C ) G (0,1) sitc/i £/ia£ gwen any / G L p (IR d ) 
wzi/i 2 — 5 < p < oo, there exists a unique weak solution to C(u) = in M. d with the property 
that u = f n.t. on M d and («)* G L p {R d ). Moreover the solution satisfies the estimate 
\\{ u )*\\p < C p H/llp, where C p depends only on d, p, \i and the constant C in ( |^.i[ ). 



The key step in the proof of Theorem 14.11 is the following lemma. 
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Lemma 4.2. Under the same assumption on £ as in Theorem 4-1. estimate ( f^.i| ) and hence 
M-8\) hold forr> 1 . 

Proof. We will prove this lemma under the additional assumption that the coefficient matrix 
A(X) is sufficiently smooth. The assumption allows us to use estimate (13.71) . For the sake 
of completeness and future reference we shall provide a version of Dahlberg's original proof 
of Lemma 14.21 in the Appendix. We point out that our proofs of the theorems stated in the 
Introduction, given in Sections 10 and 11, involve an approximation argument and do not 
rely on Dahlberg's proof. 

The case 1 < r < 10 follows easily by covering B(xo,r) with balls of radius 1. If r > 10, 
we also cover B(xo,r) with a sequence of balls {B(xk, 1)} of radius 1 such that 

B(x ,r) C \jB(x k ,l) C B(x ,r + 1) and XB(x k ,i) < C. 

k k 

Let / denote the left-hand side of (14.11) . It follows from the assumption (14.11) (with r = 1) 
and boundary Harnack inequality (12. 4p that 



/<cW \u(X)\ 2 dX <cW \u(X)\ 2 dX 

k JT(x k ,2) k JT(x k ,l) 

r r 1 r 



k 



<} [ \Vu{X)\ 2 dX < C [ [ \Vu(x,i)\- drill. 

V JT{x k ,l) JO JB(x ,r+l) 



(4.4) 



In view of Lemma 13.41 we obtain 



!<—[[ |Vd 2 dxdt < — i I \ii\-tl.rdt 



r -I -lT(x ,3.5r) ' r3 J JT(x ,3.6r) 



< —z 1 1 \u\ 2 dxdt, 



>T(x,2r) 

where we also used Cacciopoli's inequality as well as the Harnack inequality (12.41) . □ 

For 1 < p < oo, we say the L p Dirichlet problem (D) p for C(u) = in R d + +1 is solvable 
if for any / G C c (lR d ), the solution to the classical Dirichlet problem for C(u) = in 
with boundary data / satisfies the estimate ||(w)*|| p < 



Lemma 4.3. Under the same assumption on £ as in Theorem \4-l\ the LP Dirichlet problem 
for C(u) = in is solvable for 2 — 5 < p < oo. Moreover, we have ||(w)*|| p < C||/||p 
for 2 — 5 < p < oo, where 5 = 6(d, /i, Co) G (0, 1) and C = C(d, /i, C ,p) > 0. 

Proof. By the self-improving property of the reverse Holder inequality, there exists qo > 2 
such that for B = B(x , r) and Z = (z, z d+1 ) E \ T(x , 4r), 

p I \K(x, ZT < ^ K(x, Z)dx = C^l< ^±f, (4.5) 

where we used (12. 3 h and (12. 2L Since 

u(Z) = [ K{x,Z)f{x)dx, 
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one may deduce from (14.51) that 

(uy(x)<C{M(\fn(x)} 1/P0 

for any x E R , where po = q' < 2 and M denotes the Hardy-Littlewood maximal operator 
on R d . This gives < C\\f\\ p for any p > p = 2 - 5. □ 

We need the following Cacciopoli-type lemma in the proof of uniqueness. 

Lemma 4.4. Let u, v be two weak solutions in B(xo,2R) x (r/2,3r) where R>r. Suppose 
that either u or v is nonnegative. Then 

| Vw| \v\ dxdt < — / / \u\ \v\ dxdt, 

B(x ,R) r </§ JB(x ,2R) 

where C depends only on d and \i. 

Proof. Assume that u > 0. Let X E B(x ,R) x (r, 2r). Using Cauchy inequality and 
Cacciopoli's inequality, we obtain 

/ \Vu\\v\dY<-(f \u\ 2 dy) ( [ \v\ 2 dY^ 

Jb{X,D r \Jb(X,0.5t) J \JB(X,0.5r) 

c f , x , c f , M . M , 

< — sup u I \v\dY < — / |ix||v|ar, 



r B(X,0.5r) JB(X,0.5r) r J B {Xfl.br) 

where we used the Harnack's inequality in the last step. The desired estimate now follows 
by covering B(xq, R) x (r, 2r) with balls of radius r/4. The proof for the case that v > is 
similar. □ 



The next theorem concerns the uniqueness in the LP Dirichlet problem in 



pd+l 



Theorem 4.5. Let C = —div(AV) with coefficients satisfying U. ip -/ TOj) . Suppose that 
dcu z jdx = K(-, Z) E L q (R d ) for all Z E W^ 1 and some q > 1. Let u be a weak solution to 
C(u) =0 in such that (u)* E L p (R d ), where p = q' . Assume that u(x, t) -> as t -> 0+ 
for a.e. x E R d . Then u = in K+ +1 . 

Proof. Fix Z E Lg t G(X) = G{X, Z) be the Green's function with pole at Z. Choose 

<p(x) E C^(B(0,£/2)) and ip(t) E C °°(l/(2£), 21) such that <p = 1 in B{0,£/A), \V<p\ < C/£ 
and V = 1 on (!/£,£), |V^| < Ctf on (l/(2£), l/£), |V^| < C/£ on (£,2£). Then 



3 

R d+i cte,- i K d+i Sac,- 
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This gives 

i 

\u(Z)\ < C£ J' J ^ (|VG||m| + \G\\Vu\)dY 



\y\< 2 
C ^ 



t J\y\<i 



C 



1 ll <\y\< 



(\VG\\u\ + |G||Vtt|)dY 
(|VG||«| + \G\\Vu\)dY 



4 ^IWI^2 



(4.6) 



<C£ 2 F [ \G\\u\dY+% [ [ \G\ \u\dY 
Jh J \y\<* e M J\y\<t 



4/7 igih** 

* Jo J£<\y\<t 



t 



= h + I 2 + I 3 , 

where we used Lemma 14.41 for the second inequality. 

To estimate I 2 , we note that \G(Y)\ < C(Z)/£ d - 1+a for Y G B(0,£) x (£/2,3£), where 
a > 0. This leads to 

h<£zf (uydy<Cr a -i\\(uy\\ p ^0 as£^oo. (4.7) 

*- JB{Q,U) 



Next we observe that 

ii<C j 



h<c[ M(K(-, Z))M 2/e (u)dy, (4.8) 
-/|a/|<-« 

where M denotes the Hardy-Littlewood maximal operator on M. d and 

M r (u)(x) = sup {u(x,t) : 0<t<r}. (4.9) 



Since K(-,Z) G Z/(M d ), we have Ji < C HA^C^OIIp- By the assumption, 7W 2 /K U )( :E ) - > 
for a.e. x G lR d as £ — > oo. Since A^ 2 /K M ) < ( u )* e £ p (K d ), we obtain ||-M 2 /.e( M )llp ^ as 
£^ oo. 

Finally, note that as £ — > oo, 



J 3 <c/ M(K(.,Z))(u)*dy<c{ f \(uT\ p dy\ - 0. 

Thus we have proved that I\ + I 2 + ^3 — > as £ — > oo Hence w(Z) = for any Z G □ 

Theorem 14.11 is a direct consequence of Lemma 14.31 and the following theorem. 

Theorem 4.6. Let C = —div(AV) with coefficients satisfying U.1\) - I[J~^) . Suppose that the 
L p Dirichlet problem for C(u) = in R+ +1 is solvable for some 1 < p < oo. Then for any 
f G L p (R d ), there exists a unique weak solution in R+ +1 such that (u)* G L p (M. d ) and u = f 
n.t. on M. d . Moreover the solution satisfies the estimate ||(w)*|| p < C||/|| p . 
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Proof. Let / G C c (R d ) and u(X) = J Rd fduj x be the solution to the classical Dirichlet 
problem with data /. Since 

\u(x,t)\ < % [ \( u )*\dy < Cri\\(u)% < Ct-i\\f\\ p , 

1 JB(x,t) 

it follows that uo x is absolutely continuous with respect to dx on R d and K(-, X) = duo x / dx G 
LP'(R d ). In view of Theorem 14.51 this gives the uniqueness in Theorem 14.61 

To establish the existence, we let / G L p (R d ) and choose fk G C c (R d ) so that fk — > f 
in L p (R d ). Let Uk be the solution to the classical Dirichlet problem in Mi +1 with data fk- 
It follows from Lemma S3] that ||(wfc)*llp < C||/fc||p an d \\(u£ — Uk)*\\ p < C\\fi — fk\\ P - This 
implies that Uk converges to a weak solution u, uniformly on any compact subset of R!^ 1 . 
Standard limiting arguments show that — C||/|Ip an d u = f n -^- on ^ 



5 L 2 estimates for the regularity problem 

Definition 5.1. Let 1 < p < oo. We say that the LP regularity problem (R) p for C{u) = in 
W 1 ^ 1 is solvable, if for any / G CQ(lR d ), the unique solution to the classical Dirichlet problem 
with boundary data / satisfies the estimate ||iV(Vtt)||p < C||Va;/||p. 

In this section we study the L 2 regularity problem in under the condition that the 
coefficient matrix is periodic in the t direction. As in the case of the Dirichlet problem, we 
need a local solvability condition: for any xq G M d and < r < 1, 

/ \Nr(Vu)\ 2 dx < d \ [ \V x f\ 2 dx+ - [[ \S7u\ 2 dxdt\, (5.1) 

JB(x ,r) 2 KJB(x ,2r) r J JT(x ,2r) J 

whenever u G W 1,2 {T(xq 1 4r)) is a weak solution to C(u) = in T(xo,4r) such that u G 
C(T(x ,4r)) and u — f G C 1 (B(xo,4r)) on B(x ,4r) x {0}. We will also assume that the 
coefficient matrix A(X) is sufficiently smooth. However, constants C in all estimates will 
depend only on d, \i and C\ in (15.11) . 
We begin by observing that since 

|u(x,*)-u(a;,0)| < CtN 2t (Vu)(x) (5.2) 

(see |2T], pp. 461-462), condition (15. ip implies the local condition (14.11) for the Dirichlet 
problem. It also follows from (15.11) that 



I \N A {Vu)\ 2 dx<C j I \X7 x f\ 2 dx+ [[ \Vu\ 2 dxdt\. 



(5.3) 



Furthermore, if C{u) = in u G C(R d x [0, oo)) and u = f G C%(R d ) on R d , we may 

integrate both sides of (15. 3p with respect to xq over R d to obtain 



! \N A (Vu)\ 2 dx < C \ [ \V x f\ 2 dx 





\Vu\ 2 dxdt \ 


fo Jl 


i d J 



+ / / \Vu\ 2 dxdt}. (5.4) 
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Theorem 5.2. Let C = — div(AV) with smooth coefficients satisfying conditions U.l\) . U.S\) . 
U.3\) . Also assume that £ satisfies the condition ( f5.i)j . Then the L 2 regularity problem for 
C(u) = in R^ +1 is solvable. Moreover the solution u satisfies the estimate 1 1 iV( Vit) 1 1 2 < 
C||Va;/||2, where C depends only on d, fi and the constant C\ in Ii5.1\) . 

In view of (13.111) . we start with a decay estimate. 

Lemma 5.3. Under the same assumption on C as in Theorem \5.2l we have 

R \Vu\ 2 dxdt^0 as 00, (5.5) 

J Jt{0,10R)\T(0,R) 

where u is a classical solution to C{u) = in IR+ +1 such that u = f G C c (IR d ) on M. d . 
Proof. By Cacciopoli's inequality, 



H I I \Wu\ 2 dxdt < £ / / \ ill 1 (I. nit 



T(0,10K)\T(0,R) R J JT{0,UR)\T(0,R/2) 

< C I \(uY\ 2 dx 



(5.6) 



for R sufficiently large. Since {u)* G L 2 {R d ), the right-hand side of (I5.6P goes to zero as 
R -> 00. □ 

Proof of Theorem 15.21 Let / G C^R^) and u be the solution to the classical Dirichlet 
problem with data /. It follows from Remark 12.21 and (15.41) that 



(5.7) 



/ \N(Vu)\ 2 dx < c\ [ \VJ\ 2 dx+ [ \N 4 (Vu)\ 2 dx+ [ \(Q(u))*\ 2 dx 
<c{ I \V x f\ 2 dx+ [ [ \Vu\ 2 dxdt+ [ \(Q{u))*\ 2 dx 

{JR d JO JR d JR d 



To estimate the integral of |Vw| 2 over M d x (0,6), we let R — » 00 in (13. lip . In view of 
(I5.5p . we obtain 



/ f \Vu\ 2 dxdt<C [ \V x f\ 2 dx. (5i 

JO JR d JR d 



Finally, to handle the term with Q(u) in (15. 7p . we observe that u is also the unique 
solution to the L 2 Dirichlet problem for C{u) = in with boundary data /. This 

implies that (Q(u))* G L 2 (R d ) and by Theorem HZQ 

\{Q(u))*\ 2 dx <C [ \Q(u)\ 2 dx<C [ [ \Vu\ 2 dxdt 

jR d Jo JR d ^ 

<C [ \V x f\ 2 dx, 



where we used (15.81) in the last step. The desired estimate ||iV(Vw)||2 < C||Va;/||2 now 
follows from (1577) . (P|) and □ 

We omit the proof of the next lemma and refer the reader to [21] for an analogous result 
(Theorem 3.1, pp. 461-462) in the case of the unit ball. 
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Lemma 5.4. Let C = — div(AV) be an elliptic operator satisfying U.1\) - I[J~^) . Suppose that 
£(u) = in and iV r (Vw) G L p (R d ) for some 1 < p < oo and r > 0. Then u converges 
n.t. to f on R d and VJ G L p (R d ). Moreover, let 



Vi(x, t) — — — rr / / —— dyds 



t\B(x,t)\Jt J B (x,t)dyi 

dxi 



fori = l,...,d, then V{(x,t) converges weakly to in L p (R d ) as t — > 0. In particular, 



l|V,/|| p <||JV(V«)|| p . 

The next theorem concerns the uniqueness of solutions to the L p regularity problem. 

Theorem 5.5. Let C = —div(AV) with coefficients satisfying U.1\) - I[J~^) . Let u be a weak 
solution to C(u) = in such that N(Vu) G L p (R d ) for some 1 < p < j^, where a > 
is given by / TP|) . Then u = in ifu(x,t) -> as t — > 0+ for a.e. x G R d . 

Proof. The proof begins in the same way as that of Theorem 14.51 As such, this leads to 
|w(Z)| < 1% + 12 + ^3, where Ji, J 2 and J 3 are defined in (14.61) . To estimate Ji, we use the 
estimate (12. 2p and the observation that \u(x, t)\ < CtN(Vu)(x) for a.e. x G R d . This yields 
that 



C r N(Vu)(y )dy 

h }\y\<l " ' ~ " Wr" 0) — ^| 



h = C* I I \G\\u\dY<^l ^\r^ +a <Cr"\\N(Vu)\\ p 



as i — > 00, where we used the assumption p < Similarly, 

I * = T2 / \ G W u \ d V dt <M^ N{Vu)dy<Ct l - a -v\\N{Vu)\\ v 
1 J\ J\v\<i 1 J\v\<t 

which also tends to zero as i — > 00, since p < Finally, note that 

'. = 7/*/, \G\M d -f < pgs / WV„)|^0, 

as £ — > 00. Thus we have proved that Ii + I2 + ^3 — > as £ — > 00. Hence w(.Z') = for any 
Z G □ 

Recall that / G W 1,p (R d ) if / G W^?(R rf ) and V*/ G L p (M d ). The following theorem 
provides the existence of solutions with boundary data in W 1 ' p (R d ). 

Theorem 5.6. Let C = —div(AV) with coefficients satisfying U.1\) - [J~^) . Suppose that the 
LP regularity problem for C(u) = in R^ +1 is solvable. Then for any f G W 1 ' p (R d ) 7 there 
exists a solution u to C{u) = in IR+ +1 such that ||A^(Vu)|| p < C||/|| p and u = f n.t. on R d . 

Proof. Let / G W 1 ' p (R d ) and choose f k G C*{R d ) such that V x f k -> VJ in L p (IR d ). Let 
ffe be the solution to the classical Dirichlet problem in IR^ 1 with boundary data By 
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assumption, we have || JV(Vv fc ) ||„ < C\\V x f k \\ p and \\N(Vv k - S/v e )\\ p < C\\V x (f k - f e )\\ p . It 
follows that for any m > 1, V(t>fc — Vg) — » in L 2 (Q m ) as k, £ — * oo, where 

= I (z, *) e : M < 3m and — < i < 3m 
[ 3m 

Thus there exists u m G iy 1,2 (fi m ) such that v k — a(k,m) — > u m in W /1,2 (f2 m ) as k — > oo, 
where a(k,m) is the average of ffc over fi m . Clearly, C{u m ) = in fi m . By a limiting 
argument, we also have 

||Af m (Vu m )|| L p (B( o, m )) < C ||V X ./|| P , (5.10) 

where the nontangential maximal function N m (Vu) is defined in a manner similar to N m , 
but the variable t in ( 12. 9ft is restricted to (1/m) < t < m. 

Next, since w m+ i — w m is constant in fi m , one may define a function w G (M^ 4 " 1 ) such 
that u — w m is constant in VL m for any m > 1. It follows that m is a weak solution in 
Moreover, in view of (15.101) . we have 

\\N m {Vu)\\ L p m , m )) < C \\V x f\\ p . (5.11) 

Letting m — > oo in ( 15.111) gives ||iV(Vw)||p < C||Va;/||p. Similar argument also gives 

\\N(Vu - Vv k )\\ p < C\\V x (f - f k )\\ p . (5.12) 

Finally we note that by Lemma I5.4[ u converges n. t. to h G W^ p {R d ) and the average of 
over B(x, t) x (t/2, 3t/2) converges weakly to in L p (R d ) as t -> 0. This, together with 

the weak convergence of the average of |^ to |^ and estimate (15.121) . implies that jjj^ — 
on ~R d for i — 1, . . . ,d. It follows that / — fa is constant. Thus, by subtracting a constant, 
we obtain u — f n.t. on M. d . □ 

We conclude this section with two more theorems on the consequences of the solvability 
of the regularity problem. 

Theorem 5.7. Let £ = — div(A'V) with coefficients satisfying U. ij) -/ TOj) . Let 1 < p, q < oo 

arid ~ + ~ = 1- Then for C{u) = in the solvability of the L p regularity problem 

implies the solvability of the L q Dirichlet problem. 

Theorem 5.8. Let £ = — div(A^) with coefficients satisfying U.1\) - I[J~^) . Suppose that the 
L 2 regularity problem for C(u) = in IR+ +1 is solvable. Let u G W /1 ' 2 (T(x , 8r)) be a weak 
solution in T(xo,8r) such that N r (Vu) G L 2 (B(xq, 4r)) and u = f on B(xo,4r). Then 



L 



Nr(Vu)\ 2 dx<C [ \V x f\ 2 dx+- [ \Vu\ 2 dX. (5.13) 

JB(xn,2r\ r JT(xn,2r) 



B(x ,r) JB(x ,2r) 1 JT(x ,2r) 

We omit the proof of both theorems and refer the reader to Theorems 5.4 and 5.19 in 
[2T] , where the analogous results were proved in the case of the unit ball. We point out that 
Theorem 5.19 in [21] was stated for solutions in the whole domain. However an inspection 
of its proof, which extends readily to the case of shows that the conclusion holds for 

local solutions. Theorem 15.81 shows that the local solvability condition (15. ip is necessary. 
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6 LP estimates for the regularity problem 



In this section we study the L p regularity problem in for 1 < p < 2 + 5. 

Theorem 6.1. Let £ = — div(A'Vu) with coefficients satisfying U.1\) - I[TJ$) . Suppose that 
(R) 2 for C(u) = in M^ 1 is solvable. Then (R) p for C(u) = in ~R d + l is solvable for 
l<p<2 + 5. 

In the case of the unit ball, Kenig and Pipher [21] proved that the solvability of (R) q 
implies the solvability of (R) p for 1 < p < q + 5. Although it is possible to extend their 
results to the case of the upper half-space and thus obtain the LP solvability for 1 < p < 2 + 5 
from the L 2 solvability, we will provide a more direct proof. 

For the range 1 < p < 2, we follow the approach used in [13] (also see [21]) by proving 
||iV(Vw)||i < C for solutions of the L? regularity problem with LL\ t data. Our approach to 
the range 2 < p < 2 + 5, which is based on a real variable argument developed by Shen 
[25| [26] . is different from that used in [T3], [21]. The rather general approach reduces the 
problem to certain weak reverse Holder estimates. In particular, our argument shows that 
if (R) 2 for C{u) = in R d + +1 is solvable, then the solvability of (D) q for some q < 2 implies 
the solvability of (R) p , where ^ + ^ = 1. 

We start with a comparison principle. 

Lemma 6.2. Let u, v £ C(T(xo,3r)) be two weak solutions in T(xo,3r) such that u(x,0) = 
v(x, 0) = on B(xo,3r). Suppose that v is nonnegative in T(xo,3r). Then 

for any (x,t) € T{xq,t), where C depends only on d and \i. 

Proof. See Lemma 2.5 in [27]. □ 

The following is a localization result similar to Theorem 15.81 Note that here we assume 
the solvability of (D) 2 instead of the solvability of (R)2- 

Theorem 6.3. Let C be an elliptic operator with coefficients satisfying ( fi.i)) -( TOj) . Assume 
that the L 2 Dirichlet problem for C(u) = in is solvable. Then 

[ \N r {Vu)\ 2 dx<—(! \u(x,t)\ 2 dxdt, (6.1) 

where u G C(T(xo,6r)) is a weak solution of C(u) = in T(xo,6r) and u(x,0) = on 
B(x ,6r). 
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Proof. Let v(x,t) be the Green's function for C on W^ 1 with pole at X = (x Q , lOr). Fix 
x G B(xq, r ) and < t < |. By Cacciopoli's inequality, 



/ / |Vw| 2 ^g < / / |w(?/,s)| 2 cfyds 



*|-B(a;,t)| A ./ B(a . 



< 



(a:,*) fc l-^V-^) l )\ J I JB{x,2t) 

C f 2t f v(y,8) 2 , , 1 



t 3 \B(x ,t)\ Jt J B(x 



(x,2t) 



v(x ,r) 



d y ds ' ~d+r / / \ u (y^ s )\ d v ds 

r J JT(xoAr) 



< Cr ^- 2 ('!!M) // \u{y,s)\ 2 dyds 

r J JT(x Ar) 



where we have used Lemma [6.21 and v(xo,r) ~ r 1 d . Let w denote the £-harmonic measure 
for evaluated at (xq, lOr). Since 

v(x,t) u>(B(x,t)) „ x 

" TOT - 0M (w-*>s)w 

where M denotes the Hardy-Littlewood maximal operator on R d , we obtain 



, / 1 



\ 1/2 

iV^V^^^Cr^M^o^)^)^)-!^ // |u(s/,s)| 2 dj/ds . (6.2) 

ax \ r J JT(x ,4r) J 

Since the L? Dirichlet problem for C(u) = in is solvable, the reverse Holder inequality 
(14. 3 p for K = holds on B(x , 2r). The desired estimate now follows by the L 2 boundedness 
of M and (USD? □ 

Remark 6.4. Suppose that (D) p for C(u) = in R d + +1 is solvable for some p < 2. Then the 
reverse Holder inequality (14.51) holds for q = p' . It follows from (16.21) and (14. 5ft that 

V? r» / 1 r r \ V 2 

u{y, s)\ 2 dyds 



c 

< — 

r 


( 1 

V r d+l 


<c[ 


f 1 


r d+l 



1 /2 

[[ \Vu(y,s)\ 2 dyds] (6.3) 

J JT(xoAr) J 



1 



1/2 

<C(— d I \N(Vu)\ 2 dx 

J JB{x fir) 



A simple geometric observation shows that for x G B(xo,r) 



N(Vu)(x) < N r (Vu)(x) + ^ f \N(Vu)\dx. (6.4) 

r JB(x (h 6r) 



In view of ( 16.31) . this gives 

(h I WVu)\«dx) 1/q <c(±[ \N(Vu)\ 2 dx) 1/2 



(6.5) 



We will need this weak reverse Holder estimate to treat the case 2 < p < 2 + 5 for the 
regularity problem. 
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Theorem 6.5. Under the same conditions on £ as in Theorem \6.1\ the LP regularity 'problem 
in is solvable for 1 < p < 2. 

Proof. We follow the approach in [13] and prove that ||JV(Vm)||i < C for solutions of the 
L 2 regularity problem with H\ t data. To this end, let / be a Lipschitz function such that 
supp(/) C B(x ,r) and ||Va;/||oo < Cr~ d for some r > and x G M d . Let u be the solution 
of the L 2 regularity problem with data /. By the L 2 estimate ||iV(V«)||2 < C|| V x /||2, one has 
\\N(Vu)\\i,i(B(x ,Cr)) < C. Thus it suffices to show that if R > r and B(y , 10 R) fl B(xq, r) = 



5 B( 



\N{Vu)\dx <C (^-Y (6.6) 

lB(y ,R) KR/ 

for some a > 0. This follows from Theorem 16.31 and the pointwise estimate 

\u(X)\ < H/IU^^^cr)) < Cr l - d uj x (B(x ,r)) 

Cr a 



< CG(X, X ) < 



\X - X^-^ 01 ' 

where X = (x ,r). We omit the details. □ 

If Q is a bounded Lipschitz domain, it was proved in [27J that the solvability of (R) p for 
C{u) = in fl is equivalent to the solvability of (D) p i for C(u) = in Q, provided that (R) Po 
for C(u) = in f2 is solvable for some po > 1. We extend this result to the case £1 = 

Theorem 6.6. Let £ = —div(A'V) with coefficients satisfying ( fi.ij) -( TOj) . Suppose that (R)2 
for £{u) = in is solvable. Let p > 2 and q = p' . Then {R) p for £{u) = in ls 
solvable if and only if (D) q for £{u) = in is solvable. 

The proof of Theorem 16.61 relies on a real variable argument which may be formulated as 
follows. 

Theorem 6.7. Let F E L 2 (IR d ) and g G L p (R d ) for some 2 < p < q. Also assume that g 
has compact support. Suppose that for each ball B in M. d , there exist two functions Fb and 
Rb such that \F\ < \Fg\ + \Rb\ on 2B , and 

\Fb\ 2 dx < Ei sup 




B' 



1/2- 

|2 



\g\ dx 



B' 



where Ei,E 2 > and (3 > 2. Then F G L p {R d ) and \\F\\ P < C\\g\\ p , where C depends only 
on d, p, q, E\, E 2 and (3. 

Proof. This is a simple consequence of Theorem 3.2 in [26] . Indeed, it follows from the 
theorem that 

(/ \F\ p dx\ 1 < C\B(0,r)\r-% { f \F\ 2 dx\ ' +c\ [ \g\ p dx\ ' . (6.7) 

UB(0,r) J lis(0,Cr) J {J B(0,Cr) J 

The estimate \\F\\ P < C\\g\\ p follows by letting r — > oo in (16.71) . □ 
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Proof of Theorem 16.61 In view of Theorem 15.71 we only need to show that the solvability 
of (D) q implies the solvability of (R) p . More precisely, it suffices to show that ||iV(Vw)||p < 
C||Vj;/||p, if u is the solution of the classical Dirichlet problem with data / G C^M 4 ). This 
will be done by applying Theorem 16.71 to the functions F = N(Vu) and g = |Va;/|. 

To this end, we fix a ball B = B(xq, r) in IR d . Choose <p G Cq(B(xq, 8r)) such that (p — 1 
in B(xq, 7r) and \V(p\ < Cr~ l . Let A be the average of / over B(x , 8r). Write u — X = v + w, 
where v is the solution of the L 2 regularity problem with boundary data (/ — \)<p. 

We now let F B = N(Vv) and R B = N(Vw). Clearly, F = N(Vu) < F B + R B . By the 
L 2 regularity estimate, 

Wr\ I \ FB \ 2dx< ^Tml \N(Vu)\ 2 dx<£- [ \V x {(f-X) V )\ 2 dx 

\ Lt >\ J2B \ a \ JR d |-°l JR d 

<T^[ \V x f\ 2 dx<C sup ^ / \g\ 2 dx, 

\ t} \ JB(x ,8r) B'DB |-£> J B > 

where we have used the Poincare inequality. This gives the estimate needed for F B . To 
verify the condition on R B = N(Vw), we note that w is the solution of the L 2 regularity 
problem with data (/ — A)(l — tp). Consequently, w = on B(x , 7r) and we may use the 
weak reverse Holder estimate (16.51) . This leads to 

' \N(Vw)\ q dx\ 



I, 1 

\\2B 


■J 

\ J2B 


\R B \ q dx 


<_c\ 




j 

J B(x ,12r) 


<_c\ 


[w\ 


/ 

J B(x ,12r) 


<_c\ 






[\B( 


x ,12r)\. 



2B\ J 2B J 

1/2 

\N(Vw)\ 2 dx 



sl/2 (If <|l/2 

\N(Vu)\ 2 dx \ +C{— \N(Vv)\ 2 dx 

J l\B\ JB(x ,12r) 

.1/2 (If V /2 

\F\ 2 dx > + C sup < — — / \g\ 2 dx > , 

: ,12r) J B'DB { \B \ J B , J 

which gives the estimate needed for R B . Therefore, by Theorem 16 .71 we obtain ||iV(Vw)||p < 
C||V X /|| P for 2 < p < q. Finally, since the weak reverse Holder estimate (16.51) is self- 
improving, the argument above in fact gives ||iV(Vu)||p < C||Va;/||p for 2 < p < q + 5. This 
finishes the proof. □ 



7 A Neumann function on 



)d+l 

Throughout this section we will assume that £ = — div(AV) with coefficients satisfying 
(0)-(0). 



Definition 7.1. Let g G Lj oc (R d ). We call u G H// D ' C 2 (M+ +1 ) a weak solution to the Neumann 
problem, 

on 

C(u) = in Q = R d + +1 and — = g on dQ, (7.1) 

av 
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if |Vm| G L 1 (T(0, R)) for any R > 1 and 



du dip 



for any ip G Cq 



To construct weak solutions to (17. ip we introduce a Neumann function. Let = (e^) be 
the (d + 1) x (d + 1) diagonal matrix with en = • • ■ = = 1 and e^+^d+i) = —1, and 

~ ~ x f A(x,t) if * > 0, . . 

A = y4(x,t) = < , ' 7.3 
v ; \EA(x,-t)E ift<0, v ; 

where A(x,t) is the coefficient matrix for C. Let T = T(X,Y) denote the fundamental 
solution for the operator L = — div(AV) on with pole at X. Using the fact that 

L(v) = if v(x, t) = u(x, —t) and L(u) = 0, one may show that 

f (X*, Y*) = f (X, Y) for any X, Y G R d+ \ (7.4) 

where X* = (x, —t) for X = (x,t). We now define 

X(X,F) = f(X,F) + f(X,F*) for X,Y G (7.5) 

Clearly X(X, F) = X(F, X). Since \X - Y\ < \X - Y*\ for X,Y G K^ +1 , we also have 

|X(X,F)|< ( ' 



\N(X,Y)-N{Z,Y)\< 



X-Y\ d ~ v . . 

c\x-z\° (7 ' 6) 



|X — Y\ d ~ 1+a ' 



for any X, F, Z G IR+ +1 , where C > 0, a > depend only on d and \i. 

Let N(X,y) = N(X,(y,0)) for y G lR d . The next lemma shows that N(X,Y) is a 
Neumann function for the operator C on 

Lemma 7.2. faj Forl<p<2, let g G L p (M d ) and 




N(X,y)g(y)dy for d > 3 or d = 2, 1 < p < 2, 
"(A') : < , " Ka (7.7) 

(X(X, y) - X(X , y))g(y) dy ford = 2 and p = 2, 

where X = (0, 1) G T/jen u is a weak solution to the Neumann problem on M^ +1 with 

data g. 

(b) Let F G C c (R d+1 ) and 



w{X)= [ N(X,Y)F(Y)dY. (7 \t 
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Then 

I a ^^ dX = I F( P dX (7.9) 
J R i+i dxj dxi J R d+i 

for any <p G C%(R d+1 ). 

(c) If g G L 2 (M. d ) and has compact support, then the solution u defined by ( [7. 7| ) belongs 
to W 1 ' 2 (T(0, R)) for any R>1. 

Proof. Part (a) follows from 

^' ) = 2 i r a '4 {F(A ' y)} ^ AY (7 ' 10) 

for Y = (y, 0) and if G CQ(IR d+1 ). To see this, one uses 

<p(Y)= [ ^A{r(X,F)}|^rfX, 
J Rd+ i dxj 1 J dxi 

where cp is the even reflection of ip from to R d+1 , i.e., <p(X) = <p(X) if x d+l > and 
<p(X) = <p(X*) if x d+1 < 0. 
To see (b), we write 

w(X) = [ F(X,Y)F(Y)dY 



where F is the even reflection of F and deduce (17.91) from 

dw dip 



, a ij — —dX= / FipdX. 
/ K d+i oxj OXi Jm.d+1 

Part (c) follows by duality from part (b). Let h G Cq(M+ +1 ). Note that 

dv f r 8T 

— h(X)dX = 2 / g(y)dy / —(X,(y,0))h(X)dX 
d+i dxj J Rd J R d+i dxj 

= -2 / g{y)vj{y) dy, 



where 



Using 



v 3 {y)= [ T(X,(y,0))^-dX. 
J R d+i dx j 



lkilU 2 (B(o,/i)) < <?rIK'||w' 1 . 2 (T(0,.R)) < C j r||^IIl 2 (M^ +1 )' 
we may conclude by duality that | Vw| G L 2 (T(0, R)) for any i? > 1. □ 

Definition 7.3. Let 1 < p < oo. We say the LP Neumann problem (N) p for is solvable 
if for any g G L™(R d ), the weak solution u{X) = J Rd N(X,y)f(y)dy satisfies ||JV(Vu)|| p < 
C\\g\\ P - 
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The solvability of (N) p yields the existence of weak solutions of the Neumann problem 
with data in L p { 



Theorem 7.4. Let 1 < p < oo. Suppose that (N) p for C(u) = in IR+ +1 is solvable. Then 
for any g G L p (R d ), there exists a weak solution u to the Neumann problem with data g such 
that\\N(Vu)\\ p <C\\g\\ p . 

Proof. Let g G L p (M. d ). Choose g k G C c (R d ) such that g k -> g in L p (R d ). Let u k (X) = 
J Rd N(X,y)g k (y)dy. Then ||^(Vu fc )|| p < C\\g k \\ p and ||JV(Vufc - V^)|| p < C\\g k -g e \\ p . By a 
limiting argument similar to that in the proof of Theorem 15. 6[ there exists u G jy^' c 2 (lR^ +1 ) 
such that ||iV(Vu)|| p < C||g|| p and \\N(Vu k — Vu)|| p < C\\g k — g\\ p . Note that the last 
inequality implies Vu k — ► V« in L p (T(0, R)) for any i? > 1. It follows that u is a weak 
solution to the Neumann problem with data g. □ 

Remark 7.5. Suppose that (N) p for £(tt) = in R+ +1 is solvable. It follows from the proof of 
Theorem 17.41 that if 1 < p < 2, the weak solution given by (17.71) satisfies ||iV(V«)||p < C||^|| p . 

Lemma 7.6. Let u G W^(T(0,2R)). Assume that 

a *L. g-dX = (7.11) 

d+i OXj OXi 

for any ip G Cg(S(0, 2i2) x i?)). T/ien 

l»W ""Ml S ™ (^) S / m2ffl |V«r«) V2 (7.12) 

for any X, Y G T(0, where C > and /3 G (0, 1] depend only on d and \x. 

Proof. Estimate (I7.12p follows from the De Giorgi - Nash estimate by a reflection argument. 

□ 

The next theorem addresses the question of uniqueness. 

Theorem 7.7. Let u G W/ ' C 2 (M+ +1 ) and N(Vu) G L p (R d ) for some 2 < p < j^, where (3 
is given by ^7.12 ). Suppose that jjj7.11\ ) holds for any ip G C*Q(lR d+1 ). Then u is constant in 



Proof. This follows readily from the estimate ( 17. 12ft . □ 

We end this section with a result on the implication of N(Vu) G L p (R d ) on the Neumann 
data. 

Theorem 7.8. Let u G iy/ ' c 2 (IR+ +1 ) be a weak solution of C{u) = in M+ +1 . Suppose that 
N(Vu) G L p (R d ) for some 1 < p < oo. Then u is a weak solution to the Neumann problem 
on M^_ +1 with data g for some g G L p (R d ) . Furthermore, \\g\\ p < C\\N(Vu)\\ p and 

1 f 2p du 
— / a(d + i)j(x, t)— — dt — > g weakly in L p (R d ), 

P J p UXj 

as p — > + . 
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Proof. Let ip be a Lipschitz function on M. d with compact support. Define 



^ = L r a "^ iX > (7 ' 13) 

where </? is a Lipschitz function on M, d+1 with compact support such that <p(x, 0) = ip{x) for 
x G R d . Since JV(Vu) G L p (M d ), \Vu\ G D>(T(Q,R)) for any i2 > 1. Thus the integral in 
(I7.13P converges. Also note that since C(u) = in IR^ +1 , the functional A is well defined. 

Next let <f(x,t) = ip(x)rj(t) in (17.131) . where r\ is a continuous function such that r/(t) = 1 
on [0, p], r)(t) = on [2p, oo) and r\ is linear on [p, 2p]. This gives 

— / / J2a {d+1)j -^i;dxdt\, (7.14) 

which, by Holder inequality, implies that |A(^/>)| < C||iV(V'u)||p||'?/>|| J /- Hence, there exists 
g G L p (R d ) such that \\g\\ p < C ||iV(V«)|| p and 

/ a ij^~~ ~K~dX = A(ip) = f gipdx= [ gcpdx. 

JjRd+i OXj OXi J Rd J^d 

□ 



8 L 2 estimates for the Neumann problem 

To solve the L 2 Neumann problem, as in the cases of the Dirichlet and regularity problems, 
we need to impose a local solvability condition: for x G M. d , 

[ \N 1 {Vu)\ 2 dx<C 2 { [ \g\ 2 dx+ [ [ \Vu\ 2 dxdtX , (8.1) 

JB(x ,l) {.JB(x ,2) JO JB(x ,4:) J 

whenever u G W l ' 2 (T(x , 4)) satisfies with g G C{B(x ,A)) for all (p G C^{B{x ,A) x 
(-4,4)). 

The goal of this section is to prove the following. 

Theorem 8.1. Let C = — div(AV) with smooth coefficients satisfying M.l\) , U.ty) and U.3\) . 
Also assume that C satisfies conditions fl^. 1\ ) and Ii8.1\) . Then, given any g G L 2 (W d ), there 
exists a u G W l( £ (M!^ 1 ) , unique up to constants, satisfying N(Vu) G L 2 (M. d ) and \1.2\j . 
Moreover the solution u satisfies the estimate ||A^(Vm)||2 < C||y||2, where C depends only 
on d, fi and the constants Cq in ( |^.1| ) and C 2 in $8.1\) . 

In view of Theorems 17.41 and 17.71 it suffices to prove the following. 

Theorem 8.2. Under the same conditions on £ as in Theorem \8.1[ the L 2 Neumann problem 
for C(u) = in is solvable. 
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Proof. Let g G C c (M. d ) and u be given by (17. 7p . It follows from Lemma 17.21 that w G 
iy 1,2 (T(0, R)) for any i? > 1. This allows us to deduce from the condition (18. ip by an 
integration in Xq over M. d that 

/ \N A (yu)\ 2 dx < c{ [ \g\ 2 dx+ [ [ \Vu\ 2 dxdt\. (8.2) 

Note that by the pointwise estimates on N(X,Y), we have u(X) = 0(|X| 1_a! ) as |X| — > oo. 
By Cacciopoli's inequality, this implies that 



Vul 2 dxdt — > as i? — > oo. 



T(0,3H)\T(0,fl) 

In view of ( 13.61) . we obtain 

fj 

JO JR d 



Vu\ 2 dxdt <C \g\ 2 dx. (8.3) 



By dE2D, this leads to ||iV 4 (Vw) || 2 < C||#|| 2 . 

Next, using the L? bound of N^Vu) above and Lemma [5.41 we see that u — > / n.t. on 
R d and \\VJ\\ 2 < C?||iV 4 (Vw) || 2 < C\\g\\ 2 . Thus, by Remark 



\N(Vu)\ 2 dx<C\ f \V x u(^,0)\ 2 dx+ [ \N 4 (Vu)\ 2 dx + [ \(Q(u))*\ 2 dx 

I JR d JR d JR d 

<c\f \g\ 2 dx+ [ \(Q(u))*\ 2 dx\ . 

l-/R d JR d ' ) 



U) 



Finally we note that by flU}, \Q(u)(X)\ < CIX] 1 ^" for \X\ large. This implies that 
(Q(u)Y G L 2 (R d ). Since the L? Dirichlet problem for C{u) = in R^ +1 is solvable, we obtain 

\(Q(u))*\ 2 dx<C [ \Q{u)\ 2 dx<C [ [ \Vu(x,t)\ 2 dxdt<C [ \g\ 2 dx 

JR d JO J«. d JR d 

where we used ( 18. 3 p in the last step. In view of (18.41) . we have proved that [| A^(Vit) || 2 < 
C\\gh. ' □ 

We end this section with a localization theorem. 

Theorem 8.3. Let £ = —div(AV) with coefficient matrix A(X) satisfying ( fi.i)) -( TOp . Sup- 
pose that (N) 2 and (R) 2 for C(u) = in IR+ +1 is solvable. Let u G W 1,2 (T(xq, 6r)) be a weak 
solution in T(x ,6r) such that N r (Vu) G L 2 (B(x , 4r)) and |^ = g on i?(x ,4r). Then 

[ \N r (\7u)\ 2 dx <C [ \g\ 2 dx + - [ |Vw| 2 dX. (8.5) 

JB(x ,r) JB{x ,3r) r Jt{x ^t) 
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Proof. Let <p G C%°(B(x ,3r) x (-3r,3r)) such that <p = 1 in T(x ,2r) and \V<f\ < C/r. 
Let w = u — (3 where f3 is the average of u over T(x , 3r). Then for X G T(x , 2r), 

- 1. ^ « y - / <+ . y >«S ■ I dy (8 - 6) 

+ L r ail ^ {N(x ' Y)]w -^ iY 

With this representation formula, the rest of the proof similar to that of Theorem 6.10 in 
We omit the details. □ 



Remark 8.4. Let u G W^ 2 (T(x ,6r)) be a weak solution in T(x , 6r) such that N(Vu) G 
L 2 (B(x ,4r)) and g = on B(ac ,6r). Suppose that (JV) 2 and (_R) 2 for £(«) = in 
are solvable. Let E be the average of u over T(xo,4r). It follows from Theorem 18.31 that 



f4 f \N r (Vu)\ 2 dx] ' <c(4~T f \Vu\ 2 dx] 1 
WJbm J ~ \r d+1 J T ( X0 ,3r) J 

<-(-Lt [ \u-E\ 2 dx) [ \u-E\dX 

<^k[ |V«|dX<£/ \N(Vu)\dx, 

r JT(x ,3.2r) r JB(x ,4r) 

where we have used Cacciopoli's inequality, boundary regularity of weak solutions, and 
Poincare inequality. This yields the weak reverse Holder inequality, 



1/2 



C 



n , . \N(Vu)\ 2 dx\ <^/ \N(Vu)\dx. (8.7) 

By the self-improving property of the weak reverse Holder inequality, there exists p > 2, 
depending only on d and the constant C in (18. 7p . such that 

If } Vp r 1 r > V2 



|iV(Vu)| p da; f <C\— \N(Vu)\ z dx 

T JB(x ,r) J l r JB(x ,2r) 

9 L p estimates for the Neumann problem 

In this section we study the L p Neumann problem in 

Theorem 9.1. Let C = —div(AV) with coefficients satisfying U.1\ )- I[J~^) . Suppose that 
(iV) 2 and (R) 2 for C(u) = in R^ 1 are both solvable. Then (N) p for C(u) = m W^ 1 is 
solvable for 1 < p < 2 + 5. 
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Proof. In the case of the unit ball (or a star-shaped Lipschitz domain), it was proved in [21] 
that the solvability of (N) Po and (-R) Po implies the solvability of (N) p for any 1 < p < p + 5. 
We follow the approach in [T3] and [21] to treat the case 1 < p < 2. For / G L 2 (R d ) with 
compact support, let T(f) = N(Vu), where u is given by (17.71) . Since (N)% is solvable, T is 
bounded on L 2 (M. d ). With Theorem 18.31 and pointwise estimates of the Neumann function at 
our disposal, one may establish ||T(/)||i < C, where / is an FL\ t atom. The LP boundedness 
of T then follows by interpolation. 

We use Theorem 16.71 to treat the range p > 2. Let / G L 2 (M. d ) with compact support 
and F = T(f) = N(Vu), where u is given by (17. 7p . For each ball B in M. d , we choose 
F B = T(fxm) an d Rb — ^X/xrAbb) ■ Then \F\ < F B + R B . By the L 2 boundedness of T, 



[ \F B \ 2 dx<C [ \f\ 2 dx. 

JB J6B 



To estimate Rb, we use the weak reverse Holder inequality (18. 8p . This gives 



- d J b \Rb\ p oIx) <C(- d j 2B \RB?dx 



- c 0~ d L |T(/) |2 dx ) 1/2 + c L i/|2 dx ) 1/2 

Thus, by Theorem 16.71 we obtain ||T(/)|| 9 < C||/|| 9 for any 2 < q < p. 

□ 

The next theorem deals with the question of uniqueness for 1 < p < 2. 

Theorem 9.2. Let C(u) = —div(AV) with coefficients satisfying ( fi .-?)) - ( TOj) . Let 1 < p < 2. 
Assume that (R) p and (N)% for C(u) = in are both solvable. Suppose that N(Vu) G 



L p (R d ) and u G W^(R d + +1 ) satisfies Then u is constant 



in 



Proof. Since N(Vu) G L p (R d ), u = f n.t. on R d for some / G W 1 ' p (R d ) by LemmaEl Thus 
u is the solution of the LP regularity problem with boundary data /. Choose fk G Co(M d ) 
such that V x fk —>■ V x / in L p (R, d ). Let Uk be the solution of the classical Dirichlet problem 
with data fk- By Theorem 17.81 Uk is also a weak solution to the Neumann problem on R^ +1 
with data gj. G L p (R d ). Since \\gu — ge\\ P < C*||A^(VMfc — Vu^)|| p — > as k, £ — > oo, we may 
deduce that — > p in L p (IR d ). 

Next we note that ||iV(Vu fc -Vu)||p < C|| V x (/ fe - /) || p . Thus Vw fc -> Vm in L 1 (T(0, i?)) 
for any _R > 1. It follows that L d g<pdx = for any ip G CQ(R d+1 ). Consequently, we obtain 
g = 0. 

Finally since ||iV(Vufc)||2 < C\\ V x fk\\2 < oo, Uk is a solution to the L 2 Neumann problem 
with data g^. By Remark [73| is given by w fc = fi k + J Rd iV(X, y)gk(y) dy for some constant 
/3fc, if d > 3. If d = 2, one needs to replace N(X,y) by N(X,y) — N(X ,y). In any case, 
because <jr& — > in L p (M d ), we have life (a;) — (3k —* uniformly on any compact subset of 
IR^ +1 . This, together with the fact that Vu^ — > V« in L 1 (T(0,i?), implies that Vm = a.e. 
in Hence u is constant in JRl +1 . □ 
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10 Proof of Theorems II. 1L 11.21 and 



1.3 



In this section we give the proof of Theorems ll.il [L2l and [T73l stated in the Introduction. Let 

D = {(x,t) : x G R d and t > ip(x)}, 

where ip : M. d — > IR is a Lipschitz function with Lipschitz constant || V^lloo- Consider the 
bi-Lipschitzian mapping Mr^ 1 given by Y = $(X) = t) = (x, t — f>(x)). Let 

$' denote the Jacobian matrix of $ and A(Y) = (%(T)) = where 
X = $- x (y) = $(y). Note that 

J D dxjdxi J R d+i dyjdyi 

where u(Y) = u(p(Y)) and <p(Y) = tp(V(Y)). This shows that £{u) = -div(AVu) = in 
D if and only if £(u) = — div(AVu) = in Mi +1 . Consequently, the L p Dirichlet problem 
for £{u) = in D is solvable if and only if the L p Dirichlet problem for £{u) = in R^ 1 is 
solvable. The same can be said for the LP regularity and Neumann problems. 

Note that $'(*(Y)) is independent of s = and A(^(Y)) = s + ^%)). It follows 
that A(y, s + 1) = s) if t + 1) = t). Also, since 

A(y, Sl ) - A(y, s 2 ) = &(X){A(y, Sl + i>(y)) - A(y, s 2 + ^(y))}(^(X)Y , 

there exist Ci,c 2 > depending only on d and HV^IIcxj such that Ci^(p) < r/(p) < c 2 r](p), 
where r/(p) is given by (11.101) and rj(p) is defined in the same manner, using A. Therefore 
the elliptic operator £ satisfies the same assumptions as those imposed on £ in Theorem 
11.11 As a result, it suffices to prove Theorems 11.11 11.21 and 11.31 in the case D = IR^ +1 . 
Next we reduce the general case to the case of smooth coefficients. 

Lemma 10.1. Let £k = —div(AkV) and £ = — div(AV) with coefficient matrices satisfying 
< TQ]) -< rOp . Suppose that A k (X) —> A(X) for a.e. X G R d+1 as k — > oo. Let u k and u be 
solutions to the classical Dirichlet problem for £ k and £ respectively in R+ +1 with boundary 
data f G C^R^). Then, as k — > oo ; u k — > u uniformly in T(0,R) for any R > 1 and 
V-Ufc — > Vm in L 2 (R d f 1 ). Consequently, 

\\{u)*\\ p < liminf ||0fc)*ll P , 
\\N(Vu)\\ p < liminf \\N{Vu k )\\ p . 

k— »oo 

/or any p > 1 . 

Proof. We may assume that / G Cg(R d+1 ). Then 

/ AV{u- f)-VgdX = - [ AVf-VgdX 
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for any g G Cq(R+ +1 ). Recall that u(X) = J Rd fdu x = 0{\X\ l - d - & ), as \X\ -> oo. By taking 
g{X) = (u(X) - f(X))cp(X/R) where cp G C^(R d+1 ), <p(X) = 1 for \X\ < 1 and <p(X) = 
for \X\ > 2, we may deduce that || Viz U^a^a+i) < C|I^/IIl 2 (m^ +1 )- Similarly, using 

/ A k V{u k -u) -VgdX = - I (A k - A)Vu ■ Vg dX, 

for any g G Cq(R+ +1 ), we obtain 

/ \V{u k -u)\ 2 dX <C ( \A k - A\ 2 \Vu\ 2 dX. 
Jr^ 1 Jr^ 1 

It follows from the dominated convergence theorem that Wu k — > Vii in L 2 (R'J +1 ). 

To see u k — > it uniformly in T(0, we first note that u k — > u in L 2 (T(0,-R)). This 
follows from the Poincare inequality 

I \u k -u\ 2 dX <CR 2 [ \V(u k -u)\ 2 dX, 

JT(0,R) JT(0,R) 

as u k — u = on R d . However, by De Giorgi -Nash estimates, the sequence {uj} is equiv- 
continuous on T(0,R). This implies that any subsequence of {u k } contains a subsequence 
which converges uniformly on T(0, R) to m. It follows that {u k } converges uniformly on 
T(0,R)tou. □ 

Theorem 10.2. Let C k = — dzf(AfcV) and £ = — div(AV) with coefficient matrices satisfying 
the same conditions as in Lemma llO. 1\ Suppose that (i?) 2 for C k (u k ) = in R+ +1 is solvable 
with uniform estimate ||iV(Vwfc)||2 < C3II V^w/d^- Then (R) p for C(u) = in R+ +1 is solvable 
for 1 < p < po, where po > 2 depends only on d, /1 and C3. Consequently, (D) q is solvable 
for p' < q < 00 . 

Proof. By Theorem 16.11 and Theorem 15. 7\ we only need to show that (R)2 for C(u) = in 
R^ +1 is solvable. To this end let u k and u be the solutions to the classical Dirichlet problem 
for C{u) = and C k {u k ) = in R+ +1 with data / G Cg(R d ), respectively. It follows from 
Lemma [1 U . 1 1 and the uniform estimate || A^(Vw fc ) || 2 < C3 1 1 V^M^Ih that 

\\N(Vu)\\ 2 < liminf ||JV(Vtt*)|| 2 < C 3 \\V x f\\ 2 . 

k— >oo 

□ 

A similar result holds for the L 2 Neumann problem. 

Lemma 10.3. Let C k = — div(A k X?) and C = — div(AV) with coefficient matrices satisfying 
the same conditions as in Lemma \10.1\ Let 



u k (X)= / N k (X,(y,0))f(y)dy, 

JR d 

u(X)= [ N(X,(y,0))f{y)dy, 

JR d 

33 



where f G L^°(R d ) and N k (X,Y), N(X,Y) are the Neumann functions constructed in Sec- 
tion 7 for C k , C respectively. Then Vu k — > Vw in L 2 (IR+ +1 ) and consequently, 

\\N(Vu)\\ p <]immf\\N(Vu k )\\ p 

j->oo 

for any p > 1 . 

Proof. The proof is similar to that of Lemma [10. 1[ Suppose that supp(/) C 5(0, r ). 
Starting with 

AVu-VgdX = I fgdx for g e C °°(E d+1 ) 



and using u(X) = O^X] 1 d ) as \X\ — > oo, we may deduce that 



d+l 



iVufdXKCr^WfWl- 



Similarly, using 

/ A k V{u k -u)-VgdX = - I (A k - A)Vu ■ Vg dX for g e C^{R d+1 ) 

jR d+1 JR d+1 



we obtain 



./ i I ./TD>d+ 1 



V{u k -u)\ 2 dX <C \A k - A\ 2 \Vu\ 2 dX. 



Thus Vu k -»• Vm in L 2 (1R^ +1 ) by the dominated convergence theorem. □ 

Theorem 10.4. Let C k = — div(Ak¥) and C = —div(A'V) with coefficient matrices satisfying 
the same conditions as in Lemma \10. 1\ Suppose that (N)2 for Ck(u k ) = in M+ +1 is 
solvable with uniform estimate ||iV(Vttfc)||2 — C^\\du k /di/ k \\2. Then (A/") 2 for C(u) = in 
M,^ 1 is solvable. If in addition, (R)2 is solvable for C k {u k ) = with uniform estimate 
||iV(Vwfc)||2 < C4II V z Ufc||2, then both (N) p and (R) p for C(u) = in M. + +l are solvable for 
1 < p < 2 + 5, where 5 > depends only on n, /1 and C4. Furthermore, the unique solutions 
of (N) p and (R) p satisfy \\N(Vu)\\ p < C\\du/du\\ p and \\N(Vu)\\ p < C\\V x u\\ p respectively 
with constant C depending only on d, p, \i and C4. 

Proof. This follows readily from Lemma 110.31 Theorems 110.21 and 19.11 □ 

Now given a (d+l) x (d+l) matrix A(X) on R d+1 that satisfies conditions (Oil . (TP) . 
( TOD and ffTTH . we define 

A k (X)= [ A{X-Y)^ k {Y)dY 

JRd+1 

where {(p k (Y)} is a standard approximation of identity. Then A k 6 C°°(lR d+1 ) and A k (X) — > 
A(X) for a.e. X E R d+1 . Furthermore, A k (X) satisfies the same conditions ([II]), (jOl . (TOD 
and (11. lip with the same // and r/ k (p) < 77 (p), where ^(p) denotes the modulus of continuity 
for A k in the t variable. In view of Theorem 110.41 we have reduced the proofs of Theorems 
11.11 11.21 and 11.31 to the proof of the following theorem. 
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Theorem 10.5. Theorems and \1.3\ hold in the case where D = M^ +1 ; p = 2, and 

the coefficient matrix A(X) is C°° and satisfies conditions ( 17.11) . (COj), ( 11.31) and 

Proof. In view of Theorems 14 . 1 1 [5721 and IHTTj it suffices to show that the square Dini condition 
( 11. lip implies the local solvability conditions ( 15.ip and ( 18. ip . To this end, we fix xq G M. d 
and choose </?i G C^(R d ) and y? 2 G C °°(M) so that < < 1, y>i(x) = 1 in B(x ,S), 
<fi(x) = outside of B(xo,9) and < ip 2 < 1, 99 2 (^) = 1 in [—8,8], <p 2 (t) = outside of 
[—9,9]. We consider the elliptic operator L\ = — div(AiV) in the star-like Lipschitz domain 
Qo = T(xq, 20), where 

Ai{x,t) = tp 2 {t)[p 1 {x)A{x } t) + {1 - tpi{x))l] +{l-<p 2 {t))I (10.1) 

and / denotes the (d + l) x (d + 1) identity matrix. Note that A\(x, t) = A(x, t) for 
(x, t) G T(xo,8) and thus C(u) = in T(xo,r) implies Li(u) = in T(xo,r), if < r < 8. 
We claim that under the condition (11.111) . both (N) 2 and (R) 2 for L\(u) = in Qq are solvable 
and the bounding constants in the nontangential maximal function estimates depend only 
on d, /i and rj(t) in (11.101) . By the localization results for the LP Neumann and regularity 
problems in star-like Lipschitz domains (see Theorems 5.19 and 6.10 in [21]), this gives 
estimates (15. ip and (18. ip . We should point out that Theorems 5.19 and 6.10 in [21] were 
stated for solutions of L\(u) = in Qq. However the same arguments in their proofs apply 
to local solutions of Li(u) = C(u) = in T(xo, r) for < r < 8. 

To prove the claim we introduce another elliptic operator L 2 = — div(A 2 V) in Q , where 

A 2 (x, t) = <p 2 (t) [<pi(x)A(x, 0) + (1 - <pt(x))l] + (1 - ip 2 (t))I. (10.2) 

Let 

e*(X) = sn V {\A l (Y)-A 2 (Y)\ : Y e B(X,S(X)/2)}, (10.3) 
where S(X) = dist(X, <9f2 ). Since 

A x (x,t) - A 2 (x,t) = <p 1 (x) V2 (t)(A(x,t) - A(x,0)), (10.4) 

it follows that e*(X) = if \X - P\ < 1 for some P e dQ \ {(y,0) ■ \y\ < 10}. It is also 
easy to see from (TIUT4T) that if X = (x,t) with |ac| < 11 and t G (0, 1), then e*(X) < rj(2t), 
where r](-) is defined by (ll.lOp . Therefore, for P e dflo and r e (0, 1), 

KP , r) .. = Lj (f:«< c rif 4 . (1 o.5) 

r JB(p,r)nn o(X) J t 

It follows from the condition (11.111) that 

sup {h(P,r) : P G dQ ] -> as r -> 0. 

Hence, by Theorem 2.2 in [22], the solvability of (iV) 2 and (R) 2 for Li(u) = in f2 is 
equivalent to that of (iV) 2 and (R) 2 for L 2 (u) = in il . Thus it suffices to prove that the 
L 2 Neumann and regularity problems for L 2 (u) = in Q are solvable and the bounding 
constants in the nontangential maximal function estimates depend only on d and [i. We 
shall give the proof for the solvability of (iV) 2 . The proof for (R) 2 is similar. 
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Thus, let u G W 1,2 (Qq) be a weak solution to the Neumann problem for L 2 (u) = in Q 
with data g G C(dflo). We will show that 

/ \N{Vu)\ 2 da <C [ \g\ 2 da + C [ \Vu\ 2 dX (10.6) 

J B(p,i)ndn Jdn Jn 

for any P G dVL$. The desired estimate ||iV(VM)||2 < C 1 1 ^ 1 1 2 follows from (110. 6p by a simple 
argument (see the proof of Theorem 11.51 in Section 11). 

To see (110. 6p . we first note that A 2 (x, t) = I for (x, t) G flo \ T(xo,9). Hence, if P G 
dQo \ {(y, 0) : \y\ < 11}, u is harmonic in B(P, 2) fl Qo. As a result, estimate (I10.6f) follows 
readily from the solvability of the L 2 Neumann problem for Av = in Lipschitz domains 
[T8] and Theorem 6.10 in [21]. To show (110.61) for P = (y, 0) with \y\ < 11, we use the fact 
that A 2 (x,t) = A 3 (x,t) if t G [-8,8], where 

A 3 (x, t) = (p!{x)A{x, 0) + (1 - (pi(x))I. 

Let L 3 = -div(A 3 V). Then L 3 (u) = in S(x ,20) x (0,8). Since A 3 (X) is independent 
of the variable t, it follows from Theorem 8.1 in [21] (also see Theorem 1.12 in PQ) that 
both (N) 2 and (R) 2 for L 3 (v) = in R 1 ^ 1 are solvable. Consequently, estimate (110.61) for 
P = (y,0) with \y\ < 11 follows from Theorem 18.31 This concludes the proof of Theorem 
[T031 □ 



11 Uniform estimates in bounded Lipschitz domains 

In this section we give the proof of Theorems 11.41 [T751 and [L6l stated in the Introduction. Let 
fl be a bounded Lipschitz domain in IR d+1 , d > 2. Note that under the conditions (11.11) . (11.21) 
and (I1.14p . the LP Dirichlet problem for C e (u £ ) = is solvable in Q for 2 — 5 < p < 00, while 
the IP Neumann and regularity problems for C £ (u £ ) = in fl are solvable for 1 < p < 2 + 5. 
This follows directly from [12] and [21]. However, without the periodicity condition, the 
constants C (and S) in the nontangential maximal function estimates in general depend on 
the parameter e. 

The proof of Theorems 11.41 11.51 and 11.61 relies on the following two observations: 

1. The class of elliptic operators {C\ = — div(AV)} with A(X) satisfying the conditions 
in Theorem 11.41 is invariant under rotation of the coordinate system. 

2. Suppose C £ (u £ ) = in D = {(x,t) : x G R d and t > ip(x)}. Let w(X) = u £ (eX). 
Then £i(u>) = in 

D £ = {(x,t) G R d+1 : x G R d and t > tp £ {x)}, 

where ip £ (x) = e~ 1 tp(ex). 

Since HV^Hoo = HVV'Hoo, it follows from the observation (2) and Theorem II .11 that solu- 
tions of the L 2 Dirichlet problem for C £ {u £ ) = in D satisfies the estimate || (ti e )*||z, 2 (az>) < 
C|| M e||L 2 (aD) uniformly in e > 0. Similarly, by Theorems 11.21 and 11.31 solutions of the L 2 
Neumann and regularity problems for C £ (u £ ) = in D satisfy the uniform estimates 

\\N(Vu £ )\\ L 2 (dD) < C\\du £ /dv £ \\ L 2 (dD) and \\N(Vu £ )\\ L 2 {dD) < C|| V ta „M £ || L 2 (az)) , 
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respectively. 

We begin with the Dirichlet problem for C £ {u £ ) = in Q. 

Proof of Theorem 11.41 Although Theorem 11.61 implies Theorem 11.41 we give a direct 
proof. 

Let ujf denote the £ £ -harmonic measure on dfl and K £ (X, P) = dujf /da, where X G Q 
and P G dQ. Let I(P, r) = B(P, r)ndQ. We will show that for any P G dfl and < r < r , 

{-! \K £ (X,P)fda(P)X 1/2 <C[ K £ (X,P)da(P), (11.1) 

l r Jl(P ,r) ) r Jl(P ,r) 

where X G Q \ B(P , 4r) and C is independent of e. 

Fix P G dQ and < r < r . Since the class of elliptic operators {— div(AV)} with A(X) 
satisfying the conditions of Theorem 11.41 is invariant under rotation, we may assume that 

Q n B{P , CV ) = {{x, t) : x G R d and t > ^{x)} n B(P , Cr ) 

= DDB(P ,Cr ) [ " ] 

and t + to) = t) for some t > 1 and any (x,t) G IR d+1 . Without loss of generality, 
we may also assume that to — 1- Let Gn(X,Y) and Gd(X,Y) denote the Green's functions 
for L £ on Q and D respectively. Fix Y G D\ B(P ,4r). By the comparison principle (12. 5p . 
we have 

G n (X,A,(P)) G n (X,A(Po)) . _ , Pc „p > 

where A S (P) = P + (0, . . . , 0, cos). It follows that 

ug(I(P,s)) _ Gn(X,A s (P)) _ G D (Y,A S (P)) G n (X,A r (P )) 
s d ~ s ~ s ' G D (r,A r (P )) 

„ fa£(I(P,a)) G n (X,A.(Po)) 

s d G D (r,A r (p ))' 

where co>£> denotes the £ £ -harmonic measure on dD. This implies that for any P G J(Po,r), 

G n (X,A r (P )) 



K £ (X,P) <CX e , D (7,P) 



G B (F,A r (Po)) 



where K £) d(Y, P) = duj /da. Note that by Theorem 11.11 and a simple rescaling argument, 
Ke,d(Y,P) satisfies the reverse Holder inequality with exponent 2 on 7(P ,r). Thus, 



/ \K £ (X,P)\ 2 da(P) <C [ \K £ , D (Y,P)\ 2 da(P) ( 

Jl(P ,r) Jl(P ,r) V 



G n (X,A r (P )y 2 
G D (Y,A r (P )) 



< Cr- rf K(/(P ,r))] 2 f G ^X,MPo)r 2 



G D (Y,A r (P )) 



<Cr~ d [u;X(I(P ,r))} [ 
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This gives the desired reverse Holder inequality (jll.ip . □ 

Next we consider the LP Neumann problem for C £ {u £ ) = in Q. 

Proof of Theorem 11.51 By [2TJ it suffices to establish the uniform estimate ||iV(Vw e )|| 2 < 
C||g|| 2 , where u £ is the solution to the L 2 Neumann problem in Q with boundary data 
g G L 2 (dQ) and J dQ gda = 0. By dilation we may assume that \dQ\ = 1. 

Fix Pq e dfl. As in the proof of Theorem 11.41 we may assume that Q n B(Pq, CVq) is 
given by (111.21) and A(x,t + 1) = A(x,t). By Theorem 11.21 and rescaling, the L 2 Neumann 
problem for C e {u) = in D is solvable with uniform estimates. It follows from Theorem 18.31 
that 

I \N(Vu £ )\ 2 da <C [ \g\ 2 da + C [ \Vu £ \ 2 dX. 
Ji(p ,r ) Jan Jn 

By covering dfl with balls of radius tq, we obtain 

/ \N{Vu £ )\ 2 da <C [ \g\ 2 da + C I ' \Vu £ \ 2 dX. (11.3) 
Jon Jdn Jn 

To estimate f n \ Vu £ \ 2 dX, we let (3 be the average of u £ on fl. Note that 

/ \u £ - f3\ 2 da < C \ I \u £ -P\ 2 dX+ [ \ Vu £ \ 2 dx \ <C [ \Wu £ \ 2 dX, 
Jdn [Jn Jn J Jn 

where we have used the Poincare inequality. Thus 



fi / \Wu £ \ 2 da< / g(u £ - P)da <C{p) / \g\ 2 da + p / \Vu £ \ 2 dX, 
Jn Jdn Jdn Jn 

where < p < p. This implies that ||V« e ||i2(n) < Clklh- In view of All .3j) we obtain 
\\N{Vu £ )\\ 2 <C\\g\\ 2 . ' □ 

Finally we give the proof of Theorem 11.61 

Proof of Theorem 11.61 By |21j it suffices to show that 

||iV(V« e )|| 2 <C||V tari /|| 2 + i^|^||/|| 2 , 

where u £ is the solution of the L 2 regularity problem for C £ (u £ ) = in Q with data / G 
W 1)2 (dVt). By dilation we may assume that \dQ\ = 1. 

By a localization argument similar to that in the proof of Theorem 11.51 we obtain 

/ \N(Vu £ )\ 2 da<C [ \V tan f\ 2 da + C [ \Vu £ \ 2 dX. 
Jdn Jdn Jn 

Let g = The desired estimate follows from the observation 



II. 



Jn 

where we have used ||g|| 2 < C||iV(Vw e )|| 2 . □ 
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Remark 11.1. Suppose A{X) is periodic and C a (M. d+1 ) for some a > 0, Then it satisfies 
the condition (I1.14p . Furthermore, by the L°° gradient estimate in [2J, 



1/2 

2 



\vu e (x)\ < c — - / |Vu £ (r)|^r 

if = in -B(X, 2r). Consequently, the solutions of the L p Neumann and regularity 

problems in Theorems 11.51 and 11.61 satisfy the uniform estimates 

\\(Vu £ y\\ p <C\\du £ /dv E \\ p and \\(Vu E )% < C\\V tan u £ \\ p 

for 1 < p < 2 + 5. 



12 Appendix: Dahlberg's Theorem 

In this appendix we give Dahlberg's original proof of Lemma [4.21 (and thus Theorem 14.1 



Lemma 12.1. Let C = —div(AV) with coefficients satisfying U.l\) , U.S\) and U.3\) . Let u 

be a weak solution of C{u) = in Q = B(x , r) x (t — r, t + r) for some (x , t ) £ M. d+1 and 
r > 3. Then 



C ( 1 



' f \ 1/2 

\u(x ,t + l) - u(x ,t )\ < ^ J J \u(y,s)\ 2 dydsj , (12.1) 

where C depends only on d and fi. 

Proof. Let F(x, t) = Q(u)(x, t). By the periodicity assumption (11.31) . C(F) = in B(x , r) x 
(to — r, to + r — 1). It follows that 

|F(*o,*o)| < C f^I / |F(y,s)| 2 ^^ • (12.2) 

The desired estimate then follows from 

\Hy,s)\ 2 < j s \^(y,t)\ 2 dt 

and Cacciopoli's inequality. □ 

Lemma 12.2. If in addition to the assumption in Lemma \12.1\ we also assume that u = 
in 8B(xq, r) x (t — r,t + r) , then 



C ( 1 



1/2 



|«(ar,to + !)-«(*, *o)| < J [j^ J Hv, *)| 2 dyds ) (12.3) 
/or any x £ B(xq, r), where C depends only on d and \x. 
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Proof. The proof is similar to that of Lemma 112.11 In the place of (112.21) , one uses the 
boundary L°° estimate. □ 

Proof of Lemma 14.21 

The case 1 < r < 10 follows easily by covering B(xq,t) with balls of radius 1. If r > 10, 
we also cover B(x ,r) with a sequence of balls {B(x k , 1)} of radius 1 such that 

B(x ,r) C [jB(x k , 1) C B(x ,r + 1) and ^ XB{x k ,i) < C. 

k k 

Let / denote the left-hand side of (14.11) . It follows from the assumption (14.11) (with r = 1) 
and boundary Harnack inequality ( 12.41) that 



I<Cj2f \u(X)\ 2 dX < C J2H x kA) 

k JT{x k ,2) 



12.4) 



To estimate the right-hand side of (112.41) . we let v± be the Green's function for the 
operator C on the domain fix = B(xq, 2r) x (0, r) with pole at (xq, r/2). Similarly, we let v 2 
be the Green's function for C on Q 2 = B(x , 2r) x (0, 2r) with pole at (x , 3r/2). Note that 

Vi(x , r/4) > cr 2 ~ d for i = 1,2. 

Thus, by the comparison principle (I2.5p . we have 

u(x k ,l) u(x ,r/4) d _ 2 

-) 77 < C — f 77T < Cr u \ x ^ r / 4 ) 

Vi{x k ,l) Vi{x ,r/A) 

for i = 1, 2. In view of ( 112.41) . this leads to 

I < Cr 2 ^\u(x ,r/4)\ 2 Y,M x k,l)v 2 (x k ,l). (12.5) 

k 

Next we let uj denote the ^-harmonic measure on the domain evaluated at the point 
(xo,r/2). Since 



vi(x k , 1) < C / dw, 

'B(x k ,i)x{o} (12.6) 



v 2 (xk, 1) < Cv 2 (y, 1) for any y e B(x k , 1), 
the right-hand side of (112.51) is bounded by 



Cr 2 ^\u(x ,r/A)\ 2 J2 I v 2 (y,l)du. 

k JB(x k ,l)x{0} 

It follows that 

I <Cr 2 ^\u(x ,r/4)\ 2 [ v 2 (y,l)dw. (12.7) 

J B(x ,2r)x{0} 

We now consider the function 

F(x,t)=v 2 (x,t + l)-v 2 (x,t) infix. (12.8) 
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By the periodicity assumption (11. 3p . C(F) = in Q 1 . It follows that 
F(x ,r/2) = / Fdcu 



'''" J 12.!)) 



v 2 (y,l)duJ+ / F{y,r)du, 

B(x ,2r)x{0} J B(x ,2r)x{r} 

where we have used the fact that F = on dB(x ,2r) x (0,r) and F(y,0) = v 2 (y, 1) on 
B(x ,2r) x {0}. We claim that 

|F(ar 0> r/2)| < Cr 1_<i , (12.10) 
< Cr 1_d if \y-x \ < 2r. (12.11) 

Assume the claim ( I12.10I) -( I12.11I) for a moment. It then follows from ( 112.91) that 

v 2 {y,l)du<Cr 1 - d . 



B(x ,2r)x{0} 



In view of ( 112. 7p . we obtain 



C 



I<Cr d - 3 \u(x ,r/4)\ 2 <^l \u{X)\ A dX. (12.12) 

r JT{x ,t) 

Finally we note that since v 2 {x,t) < Cr 2 ~ d for (x,t) G B(x ,2r) x (0, r + 3), estimates 
(112.101) and (112.111) follow readily from Lemmas 112.11 and I12.2[ respectively. This completes 
the proof of Lemma 14.21 □ 
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